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Introduction 



Soit I'equation de Van der Pol forcee 



eil + {v? — l)u -\-u — a ^ 



ou £ est un parametre que Ton fait tendre vers 0, et ou ii designe la derivee de u 
par rapport au temps t. 

Cette equation est equivalente au systeme differentiel 



dont on a decouvert il y a plus de vingt ans certaines solutions particulieres, ap- 
pelees grands canards : pour certaines valeurs de a, dependant de £, ces solutions 

restent au voisinage de la courbe lente y — — u (obtenue en posant £ = 0) 

y compris le long de la partie instable de cette courbe, entre les deux sommets 
y — —1 ei y = \ (cf. chapitre 2). Des etudes avec les outils de I'analyse non 
standard, sur la droite reelle, ont permis de montrer que I'existence de canards 
est souvent ineluctable, mais qu'ils ont une duree de vie extremement courte : les 
valeurs du parametre a menant a des canards sont toutes dans un intervalle de 
taille exp(— M/s) (quand e — > O'*', et pour un certain M > 0). 
Ce systeme etant lui-meme equivalent a I'equation 



I'etude de ces solutions canards a pu progresser par le passage d'une variable u 
reelle a une variable dans le champ complexe [7]. La simplification de I'etude vient 
essentiellement du fait que les fonctions solutions, au lieu de n'etre que de classe 
C^, sont alors toutes holomorphes. Par ailleurs, le lien avec les points tournants 
apparait alors clairement. 

On constate d'abord que cette equation admet une solution formelle 



£u = y — 3- + 1* 
y — a — u 



dv 



(1 — u^)v + a — u , 



(E) 



00 



00 




n=0 



n=0 
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dont la propriete (exceptionnelle) est que tous les termes Vn{u) sont analytiques 
en u = 1. Ensuite, on peut chercher a faire le lien entre cette solution formelle et 
des solutions holomorphes, bornees en m = 1 quand s tend vers 0. Si on impose 
des domaines d'existence les plus grands possibles, on trouve exactement deux 
solutions v'^{u,e) et s) pour leurs parametres respectifs a'^{e) et a~{e). 

Avec ces solutions qui sont definies toutes deux sur le segment u e [—1 + p, 1], 
pour tout p > fixe, on est capable de donner une valeur a M ; on pourra choisir 

M = R{-1 + p) - R{1) = ^(^J^ - + 1)^^^^ 

Outre le premier point tournant = 1, on voit ici que I'autre point tournant u — 
— 1 joue aussi un role important. Et qu'il pourrait etre interessant de s'approcher 
de ce second point, pour remplacer —1+p par Ui = — 1 + Oe(l), Ui > — 1. 
D'ou I'idee d'etudier le voisinage des points tournants, pour montrer qu'on peut 
toujours controler la croissance, quand £ — > 0, des solutions v'^ et v~ jusqu'en un 
tel point Ui. 

Considerons le cas general d'une equation dilferentielle du premier ordre non 
lineaire 

^y' = yfix, e) + h{x, e) + ey^P{x, e, y) , 

ou e est un petit parametre complexe, ou y' designe la derivee — ^, et ou les fonc- 

dx 

tions f,hetP sont holomorphes en a; et y dans certains domaines et admettent 
un developpement en serie (eventuellement non convergent) en e. 
On peut calculer une solution formelle du type 



n=0 

pour cette equation, a I'aide d'une recurrence sur les coefficients yn{x). On peut 
ensuite au moins esperer [14] que cette serie en e, le plus souvent divergente, 
corresponde a une solution holomorphe y{x, e) quand e tend vers dans des 
secteurs axge e ]6'i, 62[ - On salt d'ailleurs qu'on a effectivement une solution, sur 
un voisinage de tout point x tels que /(x, 0) ^ 0, qui reste bornee quand e tend 
vers [12]. 

En revanche, les points Xq tels que /(xq, 0) = font partie le plus souvent [15] 
des points tournants pour cette equation : on voit bien que les coefficients y^ 
calcules sont alors en general singuliers en x — xq. l\ par ait alors difficile que 
la solution formelle y puisse etre le developpement asymptotique d'une solution 
holomorphe dans tout un voisinage de xq (la non-existence d'lme telle solution 
est precisement la definition d'un point tournant). Cependant, dans certains cas 
particuliers, il existe tout de meme une serie formelle dont tous les termes sont 
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analytiques en Xq ; ce genre de resultat est obtenii par exemple en introduisant 
un (multi-)parametre sous la forme d'lme serie formelle. Si une telle serie formelle 
existe, cela permet d'envisager I'existence d'une solution a I'equation differentielle, 
qui reste analytique en Xq. Ces solutions particulieres, holomorphes en un point 
generalement tournant pour les solutions d'une equation differentielle donnee, 
sont celles qu'on a appelees, par generalisation, solutions canards holomorphes, 
on solutions surstables. L'existence de ces canards depend souvent du choix d'un 
(multi-)parametre a{e), present dans les fonctions /, h et/ou P, qui s'ecrit lui 
aussi, bien sur, au moins comme une serie formelle en e. 

On a ainsi pu demontrer [2] dans un cas tres general que les valeurs du 
parametre a donnant des canards sont exponentiellement proches en s, et qu'on 
pent trouver un bon majorant de la constante M generalisee de la maniere suiv- 
ante : si a'^{e) et a~{s) sont deux valeurs a canards, 

\\a'^ — < C exp ( ] , pour tout £ > assez petit, 



e J 

et pour un certain C assez grand dependant de la valeur de M choisie ; comme 
dans le cas de I'equation de Van der Pol, M peut s'ecrire M = R{xi) — R{xo), 
ou R designe I'application R{x) = 3? (J^^ f{t, 0)dtj ; Xi est a nouveau un point 
dans le domaine d'existence D de plusieurs solutions canards «maximales» , point 
qu'on a interet a prendre le plus proche possible d'un autre point tournant de 
I'equation pour obtenir la meilleure inegalite. 

L'idee des recherches entreprises ici a done ete d'etudier les solutions au voisi- 
nage des points tournants pour voir si cela ne permettait pas de s'approcher de ces 
points (i.e. de prendre un xi qui tend vers une singularite quand e tend vers 0), 
et done de preciser I'inegalite ci-dessus. L'interet principal d'un tel resultat etant 
qu'il permet ensuite de donner une estimation ou au moins une bonne majoration 
des coefficients des series formelles correspondant a a. 

Dans le cas de I'equation de Van der Pol, au chapitre 2, on montrera que les 
solutions canards et existent effectivement dans des domaines contenant le 
segment [—1 + Xic^^^, 1] (au lieu de [—1 + p, 1] pour tout p strictement positif), 
ceci pour une constante Xi suffisamment grande, mais qui peut etre choisie in- 
dependante de £ ; la demonstration utilisera une equation interieure equivalente a 
(E), mais non singulierement perturbee. La propriete demontree permet ensuite 
effectivement de donner un equivalent pour les coefficients a„ des parametres a a 
canards : 

-2v^ _i/2^3^ 



n ' - nl 



" noo TTy/ne'^/^ \4 

Le resultat de I'etude des solutions au voisinage des points tournants sera 
generalise et complete au g 3.2.2 pour une equation generale que Ton reecrit sous 
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la forme 



ey' = {xPf{x) + eg{x, e, a))y + h{x, e, a) + ey'^P{x, e, a, ey) , 

dont on etudie le point tournant (pour p > 0) . 

Cette equation singulierement perturbee admet une solution (dite : exterieure) 
a distance finie du point tournant a I'interieur de certains secteurs centres en ce 
point tournant ; la demonstration de I'existence de cette fonction est rappelee 
au paragraphe g 3.2.1. Cette premiere equation est equivalente a une equation 
interieure ; celle-ci, si elle est reguliere en e, possede au moins une solution parti- 
culiere, qui pent etre connectee avec la solution exterieure. Le domaine en x de 
toute solution holomorphe bornee en e de I'equation de depart pent done etre 
etendu jusqu'a une distance Xie^^^'^^^ des points tournants : 

D = (a; e C/ |x| > Xie^^^^~^^^ et arg(a;) G ]9i, 02[} , pour des Xi assez grands . 

Dans ce nouveau domaine, cette solution n'est pas necessairement bornee quand 
£ ^ 0, mais sa croissance reste cependant controlee et ne devient pas exponen- 
tielle. 

Au g 3.3, a I'aide de ce resultat, on demontre I'existence, dans certains cas, de 
solutions surstables dans un voisinage complet d'un point tournant, si on choisit 
correctement le (p-multi-)parametre a en fonction de e. 

Au chapitre 4, on s'interesse alors a I'equation du Brusselator, transfomee en 
I'equation differentielle suivante, avec a pour parametre, 

dz{x) 2x / 1 \ a — 1 / 1\2£ 

^^^''^~d^ ^ {l + x)A^~ 2{l + xf) ~ {1 + x)^^ ~ ~ 2{1 + xfJlT^ ■ 

Le systeme du Brusselator equivalent est connu pour avoir des solutions canards 
pour certaines valeurs d'un parametre, mais I'etude du domaine de ces solutions 
dans le plan complexe n'avait pas encore ete faite. Grace a la generalisation du 
chapitre precedent, il devient facile de redemontrer I'existence des canards holo- 
morphes et de decrire leur domaine d'existence etendu. Comme pour le cas de 
I'equation de Van der Pol, on pent ensuite trouver un equivalent (quand n — > oo) 
pour les coefficients des parametres a canards. 

A noter que, pour I'equation du Brusselator comme pour I'equation de Van 
der Pol, il faut commencer par reussir a calculer des coefficients de Stokes pour 
une equation differentielle non lineaire avant de pouvoir donner les equivalents 
cites. Les multiplicateurs de Stokes, qui apparaissent autour d'un point tournant 
ne sont pas exactement calculables dans le cas non lineaire, mais dans les deux 
exemples cites, on reussit a en determiner un equivalent complet. 
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Chapitre 1 

Definitions, proprietes elementaires 



1.1 Notion de relief 

Dans tout ce qui suit, il sera souvent fait mention de nreliefz et de chemins 
descendants ce relief, notions introduites par J. Callot [3] voici quelques annees. 

Pour faire comprendre son interet, regardons I'equation differentielle singuliere- 
ment perturbee (c'est-a-dire avec e qui tend vers 0) linearisee generale 

ey' = yf{x, e) + h{x, e) + ey'^P{x, £, y) , (E.diff.) 

oil toutes les fonctions considerees sont holomorphes dans certains domaines en 
toutes leurs variables. On souhaite montrer que cette equation admet une solution 
holomorphe en ses deux variables y{x, s). Si une telle solution existe, en utilisant la 
formule de variation de la constante, on voit que necessairement y verifie I'egalite 
suivante, dans laquelle 7I designe un chemin differentiable de C allant de s a x : 

y{x, e) = y{s, ^) + ^ ^n-,^)jit,e) ^^^^ ^ ey'^P{t, e, y)) dt, 

ou F{x,e) = f{t,s)dt. Dans tous les cas envisages ci-apres, on supposera que 
F{x, s) est pen different de F{x, 0) : on aura au pire F{x, e) — F{x, 0) — 0{e) (ce 
sera bien siir le cas si / est holomorphe en £ = 0). 

Pour voir s'il existe une solution y bornee quand e tend vers 0, on s'interesse 
particulierement au terme exponentiel. Pour e reel positif, on voit bien intuitive- 
ment que si la partie reelle de F{x, e) — F(t, e) est positive pour un certain t, 
I'integrale est en general exponentiellement grande quand s tend vers 0. C'est 
a cause de cela que Ton considere I'application relief qui a tout a; de C asso- 
cie R{x,e) — ^{F{x,e)) dans R, dont les lignes de niveau out pour equation 
^{F{x,e)) = constante. II parait alors naturel de supposer I'existence de solu- 
tions dans des domaines {x/ 37^, chemin descendant le relief en chacun de ses 
points}. En fait, on en demande un pen plus : 

Definition 1 On dira d'un domaine D qu'il est accessible si : 
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1. pour tout X E D on peut trouver un chemin differentiable Jxi'^) allant d'un 
point fixe s e 2) (appele sommet deD) ax; 

2. 7^ est tel que ^^^^^ < -C^,\F'{-fx{u))i^{u)\ pour tout u ; 

3. les C-y^, tons positifs, admettent une home inferieure strictement positive : 
inf C^^ existe et est eqal d C > 0. 

Si on a une telle propriete, le terme exponentiel de I'integrale sera exponentielle- 
ment petit quand e tend vers 0. et il est alors raisonnable d'esperer pouvoir 
controler la valeur de I'integrale pour tons les a; G D. 

Si £ tend vers dans une autre direction (arg(£) = ^ 7^ 0), on change de 
relief, et done de domaine accessible, en optant cette fois pour 
condition d'accessibilite restant la meme avec le nouveau relief. On prefere en 
general effectuer les etudes avec un relief fixe (qui donne des chemins fixe) ; 
c'est la raison pour laquelle on travaillera toujours avec a.Tg{e) constant et \e\ 0. 

Pour I'integration numerique aussi, la propriete de relief est importante : il 
n'est pas difficile de verifier que pour deux solutions de I'equation differentielle 
et telles que i/aixo) 7^ ybi^o), la difference entre les deux fonctions se comportera 
de maniere totalement differente (pour e arbitrairement petit) suivant qu'on part 
de Xq le long d'un chemin qui monte, on qui descend du relief. En effet, on aura 

1 f F(x,s)-F(t,e) / \ 

{ya-yb){x,£) = {ya-yb){xo,£) + - I e ^ l£dP[t,e,ya,yb)]dt . 

L'integrale dans cette egalite sera, si ^y^ descend le relief, exponentiellement petite 
en £ ; ce qui signifie que les deux solutions resteront voisines (du moins si le chemin 
entre xo et x est de longueur finie). Mais I'integrale deviendra exponentiellement 
grande si ce chemin monte le relief ; et done dans ce cas ya et yb s'eloigneront tres 
rapidement I'une de I'autre. 

Lors d'une integration numerique d'une equation differentielle du type cite, 
en prenant un petit e, il sera done tres important de partir d'un sommet Xq, 
«haut » sur le relief, et d'en rester aux points accessibles a partir de ce point. 

Cette constatation nous amene a deux nouvelles definitions : 

Definition 2 On appellera montagne, pour le relief R{x), par rapport au point 
Xq E C (dans les cas cites, Xo sera un point col pour le relief R), un domaine de 
C connexe et simplement connexe tel que xq est accessible avec ce relief a partir 
de tout element du domaine : R{x) > R{xo). 

De meme, on appellera vallee un domaine, toujours simplement connexe, ac- 
cessible a partir de xq ; les points x du domaine verifient done : R{x) < R{xo)- 
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1.2 Serie formelle Gevrey 

Quand on essaye de donner une solution sous forme de serie a une equation 
differentielle singulierement perturbee, on arrive a des series divergentes, mais 
dont les termes verifient cependant des proprietes de croissance. On est amene 
naturellement a introduire les series Gevrey (cf. [14], [15]). 

Definition 3 On notera dans toute la suite le secteur ouvert centre en Xq, de 
rayon p et d'angle ]6i, sous la forme 

S{xo, p, Oi, 62), soit |x e C y^O < |x — Xo| < p, et arg(x — xq) e ]9i, ^2[| 

Si le rayon p n'est pas precise, il sera suppose infini. 

Par extension, on parlera aussi de secteurs centres en I'infini : 

S{oo, g, 01, 62) = ^xeC ^ \x\ > g, et arg(a;) G ]9i, ^2[} 

Definition 4 On dit d'une serie formelle 

00 

n=0 

qu'elle est Gevrey d'ordre 1 si la serie a un rayon de convergence nul, mais qu'il 
existe des constantes positives A B telles que, pour tout n 

\an\ < AB'^n! 

Cette definition s'etend bien sur aux series de fonctions 



n=0 



en remplagant |a„| par ou || • || designe la norme que Ton choisit pour les 

fonctions 

Definition 5 On dit qu'une fonction a{e) admet la serie ^ fln^" comme developpe- 
ment asymptotique Gevrey dans un certain secteur S — 5'(0, p, ^1, ^2) si a{e) 
est holomorphe sur S et s'il existe A et B, constantes positives, telles que : 



pour tout N assez grand et tout e dans S , 



N-l 



a{s) - X ane"" 



n=Q 



Dans le cas d'une serie de fonctions (de la forme ci-dessus), si la majoration est 
independante de x, pour tout x dans un certain domaine D, on parle de serie 
asymptotique uniforme sur D. 



12 



On sait qu'il existe une infinite de fonctions holomorphes a(s) admettant 
'^QnS"' comme developpement asymptotique en (on ecrit a ~ ^a„e"). 

Si a ~ ^QnS"', on sait par ailleurs qu'on obtient une approximation exponen- 
tiellement bonne en e de la valeur de a{s) en sommant la serie des a„£" jusqu'a 
ce que, a e fixe, le terme la^s"! soit minimum (on parle alors de sommation «au 
plus petit terme » ). 

1.3 Serie formelle et vraie solution 

On se place dans le cas oil une equation differentielle (singulierement pertur- 
bee) du type (E.diff.) 

£y' = yf{x, e) + h{x, e) + ey^P{x, e, y) 
admet une solution formelle 

oo 
n=0 

Definition 6 Si la serie ^2/n(a;)£" est solution formelle de (E.diff.), on ap- 
pelle dassiquement ncourbe lentez la fonction yo{x), obtenue pour la valeur 
du parametre £ = 0. 

Le probleme est de montrer qu'il existe une solution holomorphe y(x, e) (ex- 
istant dans certains secteurs centres en pour la variable e, et dans certains 
domaines pour la variable x) , qui admette effectivement y comme developpement 
asymptotique quand e tend vers 0. 

La reponse generale est loin d'etre evidente. Un certain nombre d'articles ([9], 
[6]) out ete recemment publies a ce sujet, partant de I'equation de croissance 
cristalline 

^2^/// + = cos $ 

avec les valeurs aux bornes 

$(s) — > — 7r/2 quand s — > — oo, $(s) — > 7r/2 quand s — > +oo. 

lis demontrent que s'il existe bien une solution formelle $ = J2'^=o^^^^n{s), ou 
$0 verifie bien les valeurs aux bornes et limoo = pour tout n, il n'existe 
pas de solution holomorphe a I'equation verifiant les limites voulues en I'infini (le 
probleme venant d'un terme exponentiellement petit en e, qui n'apparait done 
pas dans la serie). lis donnent aussi des ebauches de procedures pour determiner, 
pour des equations plus generales, si on est en presence de tels cas. 

Le travail principal dans cette these correspond au probleme inverse de celui 
envisage pour I'equation de croissance cristalline. Sachant qu'on a une solution 
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formelle d'une equation differentielle qui verifie une certaine propriete (dans notre 
cas, la continuite en un point qui est en general singulier), il s'agit de montrer 
qu'il exists effectivement une vraie solution a I'equation admettant cette serie 
formelle comme developpement asymptotique. Cette solution verifie alors la meme 
propriete (la continuite en ce point, ou plus exactement le fait que cette fonction 
reste bornee au voisinage de ce point quand e tend vers 0). 

Meme pour ce probleme particulier, on connait des exemples ou I'existence 
de solutions formelles a fonctions coefficients continues en a; = n'implique pas 
I'existence de vraies solutions bornees, quand e tend vers 0, au voisinage de a; = 0. 
En effet, pour I'equation 

su — xu-\- exp I — = 



la solution nulle est solution formelle de I'equation, qui est evidemment analytique 
en 0. Cependant, les vraies solutions s'ecrivent (en utilisant la formule de variation 
de la constante, qui met a nouveau le relief en evidence, fig. 1) 

u{x) = u{xo) exp (^^) + £ I exp exp dt 

ne peuvent pas etre bornees en x. independamment de e, dans un intervalle ouvert 
autour de a; = 0. Ceci etant vrai quel que soit les xq et u{xq) choisis : pour Xq > 0, 
la solution sera bornee en e pour tout a; > 0, mais ne le sera pas pour un a; < 
quelconque fixe. On remarquera cependant que, dans cet exemple, les coefficients 
de I'equation ne sont pas tons holomorphes en e. 



1.4 Points tournants d'une equation differentielle 

On trouve diverses definitions des points tournants pour une equation differen- 
tielle dans la litterature. La plus classique [14] pour une equation du type (E.diff.) 
est la definition npar defautz suivante : 

Definition 7 Xq est appele point tournant, pour une equation (E.diff.) admettant 

une solution formelle y = X]l/n^"; s'il n'existe pas de vraie solution de I'equation 
admettant comme developpement asymptotique la serie formelle y, uniformement 
pour tout x dans un voisinage de xq (quand e tend vers 0). 

On remarquera que les points tournants sont bien des points du plan complexe 
en X : les series formelles y, sauf si elles convergent, ne sont des developpements 
asymptotiques de solutions de I'equation que pour des e tendant vers dans 
certains secteurs en e centres en 0. 

D'apres ce qui a ete dit a propos du relief, on constate que les seuls points 
susceptibles d'etre des points tournants sont les points cols du relief donne par 
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Fig. 1: Lignes de relief avec R{x) = 3? [f^ tdt) ; les signes — et + designent respective- 
ment les deux vallees et montagnes par rapport au point col 0. 

^(^J^ f{t)dt^ , autrement dit la ou / s'annule. Ce qui est le cas dans I'exemple 

precedent (fig. 1). Ou bien evidemment des points singuliers de ce meme relief. 
On peut voir en efi^et que pour un point normal x, il existe toujours un point nau- 
dessusz de x a partir duquel un voisinage complet de x est accessible ; et il existe 
done une solution bornee dans ce voisinage. En revanche, pour un point Xq qui est 
un point col, ou singulier, le developpement asymptotique sera valable en general 
au mieux sur une montagne naissant en Xq et I'essentiel du domaine accessible a 
partir du sommet de cette montagne. Pas forcement tout le domaine accessible, car 
la validite du developpement jusqu'en xo lui meme ne soit pas forcement assuree. 
Dans ce qui suit, nous etudierons comment se comportent les solutions d'lme 
equation differentielle au voisinage des points col pour le relief correspondant. 
Nous verrons que certains de ces points ne seront pas, exceptionnellement, des 
points tournants suivant la definition donnee ci-dessus. 
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1.5 Phenomene de Stokes 



Autour d'un point tournant, correspondant a une solution formelle unique, 
on trouve en general differentes solutions existant dans des (sous-ensemble de) 
secteurs ouverts centres en ce point ; on retrouvera cette propriete classique au 
theoreme 10. Ces domaines d'existence sont tous differents, mais I'intersection de 
deux de ces domaines distincts n'est pas forcement vide. L'explication de cette 
difference de comportement de deux solutions, tres voisines pour certains x, quand 
on s'approche des limites d'un domaine (I'une des solutions reste bornee, alors que 
I'autre y devient singuliere) est le phenomene de Stokes. 

Dans le cas lineaire, oii les domaines d'existence sont exactement des secteurs 
ouverts, on salt [14] que la difference entre deux solutions est un terme du genre 
exp {—P{x)/e). Ce terme est exponentiellement petit en e dans I'intersection des 
domaines d'existence (ce qui explique que le developpement asymptotique reste 
le meme), mais il devient exponentiellement grand quand x atteint les limites des 
secteurs : pour £ > 0, quand la partie reelle de P{x) devient negative. 

Dans le cas non lineaire, le phenomene est certes plus complexe, mais reste 
analogue. Quand on parlera ici de calculer un coefficient de Stokes, il s'agira de 
determiner la difference de deux solutions a une equation donnee, au moins pour 
certains x oil ces deux solutions sont toutes deux definies. Cependant alors que 
dans le cas lineaire, on pent trouver le terme exact, on n'aura dans le cas general 
qu'un equivalent (dans le meilleur des cas). 

Le phenomene de Stokes est indissociable des points tournants : s'il n'est 
pas present, c'est qu'il existe une solution holomorphe univalente tout autour 
du point tournant, correspondant a une unique serie formelle. Dans ce cas, si la 
serie formelle est bien definie en ce point, la solution est holomorphe jusque la. 
On n'a plus alors de point tournant, par definition. 

1.6 Solutions ncanardsz 

C'est dans ce cas particulier ou il existe une solution holomorphe jusqu'en 
un point col du relief, qui est un hcandidat point tournantz, que I'on parlera de 
solution canard : 

Definition 8 La (vraie) solution d'une equation dijferentielle du type (E.diff.) 



dont le relief correspondant J f{t, e)dtj admet un point col est appelee solution 

canard si ce point col n'est pas un point tournant pour elle. 

On considerera, dans la suite, uniquement les canards dans des domaines ou- 
verts du plan complexe. lis seront bien entendu holomorphes dans ces domaines. 
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Sur la droite reelle, il est en revanche possible de trouver des canards moins 
reguliers, dont seules les premieres derivees restent bornees dans un intervalle 
autour du point col. 

L'existence de solutions canards a deja ete demontree dans des cas tres generaux. 
Parmi les resultats recents, on pent citer les deux theoremes suivants : 

Theoreme 1 [8] Soit V equation 

eu' = f{x)u + eP{x, u, e) 

ou les fonctions f et P veriji,ent les proprietes suivantes : 

1. La fonction f est analytique dans un voisinage complexe d'un intervalle reel 
[a, b] autour de 0, et est a valeurs reelles sur cet intervalle. 

2. On a xf(x) > pour tout x reel non nul de ce voisinage; de plus, il existe 
A > et p entier impair tels que f(x) — Xx^ (1 + 0(x)) au voisinage de 0. 

3. Le fonction P est analytique dans un voisinage de [a, b] x {0} x {0} dans 

Alors on a I'equivalence suivante : il existe une solution formelle u = 
dont tous les coefficients Un sont analytiques au voisinage de si et seulement si 
il existe une solution canard, tendant vers Uq quand e 0, uniformement dans 
un voisinage complexe de x — 0. 

Ce premier theoreme montre que dans le cas des solutions canards, il y a en 
fait une relation presque automatique entre l'existence d'une solution formelle 
et l'existence d'une vraie solution. Le cas des equations avec parametre est un 
pen different, puisqu'on ne trouve, formellement, un parametre que sous la forme 
d'une serie formelle. Dans le cas de ces equations, on connait le resultat suivant : 

Theoreme 2 [2] On considere I'equation avec le multiparametre a suivante, definie 
pour X dans un domaine T) : 

£u' = ^(x, u, a, e) , 

oil ^ est une fonction analytique dans un ouvert connexe Dq, C x C x x C. 
On suppose que 

1. II existe G C*' et une fonction analytique Uq tels que "^{x, Uo{x), Oq, 0) = 
pour tout X E T) (il faut bien sur que si x E D, (x,Uo(x),ao,0) soit dans 

2. La fonction f{x) — ^a;, mo(x), ao, Oj admet un unique zero, xq dans T), 

qui est d'ordrep. Cette fonction f est celle donnant le relief pour I'equation, 
avec la relation R{x) — 3? (^J^ f{t)dt) ; xq est done le seul point col pour ce 
relief dans T). 

3. Le domaine Ti \ {xq] est la reunion de p + 1 domaines accessibles, ces do- 
maines etant eventuellement non homes. 
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4- On note : x — > 'if{x,uo{x),a,0). L' application 
J : CP — 

a ^(vI/„(a;o),KK),...,*i^-'H^o)) 

est un diffeomorphisme local au voisinage de oq. 

Alors I'equation dijferentielle admet une solution canard u(x), voisine de Uq, pour 
une valeur du pararnetre a tendant vers avec e. 

Ce dernier theoreme d'existence sera redemontre au g 3.3 (d'une autre maniere 
que dans Particle cite). 
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Fig. 2: Canard, dit nde Van der Polz 
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Chapitre 2 

Les canards de Van der Pol 



2.1 Introduction 

L'equation de van der Pol [13] 

eii + {u^ — l)u -\-u — Q 

est un exemple connu d'equation differentielle dont les solutions presentent des 
oscillations rapides entre deux periodes de relaxation. 
On considerera ici l'equation de van der Pol «forcee» 

eu + {u^ — l)u + u = a 

ou ii designe la derivee ou £ est un parametre arbitrairement petit et a un 
parametre reel. 

Dans toute la suite, on se limitera au cas a > (le cas a negatif lui etant 
symetrique). 

Pour etudier les solutions de cette equation, on se place dans le plan des phases 
(ou plan de Lienard, [10]), en prenant le systeme equivalent suivant : 

{eii — y — ^ + u 
y — a — u 

Une etude sommaire a partir de ce systeme montre que, pour a > 1, la trajec- 
toire rejoint rapidement le voisinage de la cubique d'equation y = — u, puis 
joint le point A(a,a'^/3 — a) qui est stationnaire. Pour a <1, a pas trop proche 
de 1, la trajectoire finit par contre en un cycle comprenant deux parties stables 
de la cubique et deux segments horizontaux (cf. figure 3). Une etude locale au 
voisinage de A montre d'ailleurs qu'il y a une bifurcation de Hopf au voisinage de 
a = 1 : le point stationnaire A est attractif si a < 1, et repulsif si a > 1 et on a 
dans ce cas un cycle limite dont la taille est d'ordre \/l — a (pour les a tres pen 
differents de 1). 
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Fig. 3: Cycle pour a = 



Le phenomene dont il est question ici est cependant distinct, quoique tres 
voisin quand on fait varier a, de celui de la bifurcation. II a ete decouvert quand 
E. Benoit, J.L. Callot, F. et M. Diener ont cherche a voir, d'abord numeriquement, 
comment on passait (pour a < 1) d'un cycle de petite taille ^/l — a au cycle 
beaucoup plus grand de la figure 3 ; d'apres la propriete de dependance continue 
des solutions (et done du cycle) par rapport au parametre a (qui intervient de 
maniere non singuliere dans I'equation), on savait qu'il devait exister un regime 



21 



intermediaire entre ces deux cycles. 

3i 




-3-2-1 1 2 

X 

Fig. 4: Canard intermediaire; B, d'abscisse 5, est appele ncolz du canard 
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Fig. 5: Grand canard 
Les premiers resultats obtenus sur ce probleme sont les suivants : pour e arbi- 
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trairement petit fixe, on montre [1] que quand a varie au voisinage de 1, on pent 
obtenir un regime intermediaire entre les deux types de trajectoires cites, dans 
lequel la trajectoire suit la cubique sur un morceau (aussi grand que Ton veut) 
de la partie instable de cette courbe ; on appelle les solutions correspondantes des 
« canards » . Plus precisement, 

1. pour des valeurs a voisines de 1 bien choisies, la trajectoire suit la cubique 

entre ses deux sommets (l,— |), (— 1, |) : on a alors ce qu'on appelle des 
{(grands canards )> (cf. figure 5) ; 

2. la trajectoire suit la cubique jusqu'au point d'abscisse 6 e ] — 1, 1[ si et seule- 

ment si |a — q;| = exp[— + 0£(1))] , ou a est une des valeurs a grand 

canard, et R{h) = J^{^ — + (voir figure 4, et la figure 2 pour 

I'origine de I'appellation de « canard » ) ; 

3. si a est une valeur a grand canard, la trajectoire passe aussi dans un voisinage 

de taille s du point {u — l,y = —2/3) et la solution du systeme correspon- 
dante est done necessairement continue en ce point (ou plutot : bornee en e 
dans un voisinage de ce point). 

Ce sont ces solutions particulieres qui vont etre etudiees. Dans ce but, on 
efi^ectue le changement de variables suivant, pour se placer au voisinage de la 
cubique : y — ^ — u + ev, ce qui nous mene au systeme 

{ii — V 
ev — v{l — u^) + a — u 

On elimine la variable temporelle : 

ev^ = {1 - u^)v + a - u . (E) 

Comme courbe lente pour cette equation (i.e. pour le parametre e = 0), on 
trouve la fonction u i— > —fz:^, qui n'est continue en 1 que si a = 1. Dans ce cas, 
on trouve facilement la solution formelle 

oo oo 

Q; = ^ v(u) = ^ Vn(u)e'' 

n=0 n=0 



avec les formules de recurrence suivantes 



vo{u) 



ao = 1 Un+l = J2 ^j(l)^n-j(l) 

j=Q 



(2.1) 



u + 1 



Vn+l{u) = Vo{u)- 



u-1 
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Une etude de cette serie, realisee par M. Canalis-Durand avec les outils de 
I'analyse Gevrey, a montre qu'elle etait divergente, mais Gevrey d'ordre 1 [5]. 

Dans [7], en etudiant I'equation dans le champ complexe, il est demontre que : 

1. pour e e R"*", il existe une unique solution tel que la solution 

existe en particulier, pour n'importe quel p reel positif, dans I'ensemble du 
secteur infini S{-1 + p, -f , ^) = [z e C / arg(;2 + l-p)e]-f,^[}; 
v'^ correspond dans son domaine de definition a la solution formelle v ; 

2. pour a~ — «+, on obtient aussi une vraie solution v~ existant dans un domaine 

symetrique de celui de ; 

3. il existe une solution {a{e), v(u,e)) holomorphe en e avec une singularite en 

6 = qui prolonge la solution [q.~^,v^) dans un large secteur centre en 
(pour la variable e) ; cette solution verifie en plus la propriete : a~{e) = 
a (ee"^*''), v^{u,6) = v {u,se~'^'''^). 

4. pour tout £ e C assez petit tel que | arg(£)| < | + 5, on a 

\a{se-''^) - a{e)\ < exp(^^(^tl±fl^^ 

ou 5 et 5' sont des constantes positives, et R{x) = J^{t — l){t + lYdt. 

L'analyse a partir de cette serie formelle ne permet pas cependant de voir ce 
qui se passe au voisinage de (—1), puisque la courbe lente y a une singularite; 
alors que le dernier resultat ci-dessus laisse penser qu'on gagnerait en precision si 
on pouvait remplacer 5' par quelquechose qui tend vers avec e. 

Nous allons regarder specifiquement ce qui se passe en ce point en introduisant 
un nouveau changement de variable qui cree une loupe pres de (—1) : 



u 



d'ou 



dY £^/^ 
(E) <^ Y— = 2XY{1- —X) +a+l-e^/^X (E') 
dX 2 

La courbe lente (£ ^ =r- a ^ 1) suivra I'equation 

dY 

Yo-^ = 2XYo + 2 (2.2) 

Nous allons etudier les solutions de (2.2) bornees a I'infini et montrer qu'il ex- 
iste une vraie solution de I'equation differentielle (E') qui admette un developpe- 
ment asymptotique du type Y{X, e) = X^^o Yn{X)e^^^ . Cette solution correspon- 
dra done aussi a la serie formelle obtenue de (2.1) par le changement de variable ; 
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I'intersection de son domaine d'existence avec le domaine de z'^ etant non nulle 
et comme elles sont egales dans cette intersection, par unicite, Y{X,e) sera le 
prolongement analytique de z'^ jusqu'au point —1. 

L'existence de cette solution permettra de trouver un equivalent exact de la 
difference a'^ — a' (paragraphe 2.4). A I'aide de cet equivalent, nous en deduirons 
ensuite un equivalent pour les coefficients a„ de a (dans le paragraphe 2.5). 

2.2 Etude de la courbe lente 

II est possible de donner, d'une certaine maniere, une solution exacte de I'equa- 
tion de la courbe lente (2.2) qui verifie les conditions aux limites voulues. Plus 
exactement, nous allons demontrer le theoreme suivant : 

Theoreme 3 II existe des solutions de V equation differentielle 

rIY 

(2.2) Yo-^ = 2XYo + 2 

qui tendent vers en I'infini et dont les domaines d'existence contiennent des 
secteurs ouverts d'ouverture au moins ^ en I'infini. En particulier, il y en a 

une qui est holomorphe dans un secteur S ( oo, g, , — ] , et une autre dans un 

\ 6 6/ 

oo, ^, ,— ). Ces deux solutions peuvent 

6 6J 

etre prolongees jusqu'en 0. 

Pour demontrer ce theoreme, on commence par effectuer le changement de 
variable suivant : 

yo(^) =X^ + z{X) =^ (2.2) z'{X) 



X^ + z{X) 

On regarde ensuite la fonction reciproque de la fonction X i— > z{X) : on considere 
en fait z i— > X{z) ; on salt, d'apres les theoremes sur les fonctions implicites que 
si z' ne s'annule pas, la derivee de X s'ecrira : 

dX(z) _ X\z)+z 
dz ~ 2 ■ ^^-"^^ 



On pose ensuite 



u' 1 
X — —2 — , d'ou u" + -uz — 0, qui est une equation d'Airy. 

u 4 

Les solutions de I'equation d'Airy peuvent toutes s'ecrire (par exemple) a partir 
des fonctions d'Airy Ai et Bi. 
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On note 

Une des solutions correspondra par exemple ku = A\{—iJ,pz 

A\[-lifz) 

Or, les resultats connus (voir [11]) sur Ai sont : 



(2.4) 



oil ^ = f-^^^^ i^^ prenant la valeur principale de la racine), ceci pour 

— TT + 6 < arg z <% — S {S > 0) . 

Comme arg(— /xj^z) = arg(2;) + 7r/3, on s'interesse aux z tels que — ^ + S < 
argz < ^ + S, et on voit que dans ce cas 

^-z + O 



Cette derniere egalite permet d'affirmer, en utilisant le theoreme d'inversion lo- 
cale, que la fonction z ^ X{z) a un inverse quand z est dans un voisinage de 
I'infini dans le secteur propose, soit la fonction X i-^ z{X). 

Et alors Yo{X) ^ X^ + z{X) = 0{l/^/\z\) tend vers quand X ^ oo pour 

-l + S'<argX<^-S' 
6 6 

On obtient bien une solution Yq bornee sur un large secteur de sommet O. 
D'apres I'equation difFerentielle (2.2), Yq ne pourrait tendre vers I'infini que quand 
X tend vers I'infini. En effet, supposons qu'il existe une suite de points {Xn) 
convergeant vers un point fini X^, telle que Yo{Xn) soit une suite non bornee. 
Alors la fonction uq — I/Yq verifie I'equation 

^ = -2Xul - 2ul , 
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et elle est telle que mo(X„) 0. Comme uq n'a pas de singularite en Xqo, et que 
d'apres le theoreme de Painleve elle verifie mo(-^oo) = 0, 11 s'agit de la I'unique 
solution de son equation qui s'annule en Xqo, done de la fonction nuUe. Ce qui est 
absurde. Les seules singularites de Iq sont done des zeros de Yq. 
L'equation differentielle nous indique aussi que 

comme Yq est bornee dans un domaine comprenant un secteur en I'infini, d'apres 
la formule de Cauchy, Y^ est bornee elle aussi a I'interieur de ce domaine. Done 
quand X devient grand, on a necessairement que 2XYo + 2 0. D'ou le resultat. 

Si on veut aller plus loin dans le developpement de Yq, il sufiit de decomposer 
la fonction Yq = —1/X + Y. La fonction Y est solution (bornee quand X devient 
grand dans le bon secteur) de 

~ 1 Y Y' 

Y'Y = 2XY + — ^ + — . 

X3 X^ X 

Y = o{l/X) ; avec Cauchy, on en deduit done que Y' est au plus un o(l/X^). 
Dans l'equation differentielle, les deux termes les plus importants sont cette fois 
2XY et l/X^, done Y - -1/2X^ 
Au total, 

^ _ 1 1 f 1 
^'-~X~2X~^^''[x~^ 

On pent montrer que la solution proposee existe aussi sur I'ensemble du demi- 
axe {z E C / argz = — f }■ 



En effet, pour z — —jt (t E 

X{z) = 2//j2_^ argX(^) = argj^ + tt - = - 
dX dXdz X'^(z)+z, jV. oAi'^O^^) 



dt dz dt 2 ^ 2 V Ar(/xt) 



Le calcul de la derivee de 4/i^Ai'^(/xt) — tA\^{nt), montre qu'elle est nulle ; cette 
fonction est par consequent une fonction constante, et reelle, done -jt^ est une 
fonction reelle qui garde un signe constant (puisqu'on I'obtient de la precedente 
en divisant par Ai^(//i) e R+). 

D'ou Ton pent deduire que t i— > ^ est reelle et uniforme sur W^, done in- 
versible. 
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On peut bien sur prolonger analytiquement cette fonction sur tout domaine 
simplement connexe ne contenant pas de singiilarite de Yq (les seules singularites 
possibles etant les points ou Yq s'annule). Cela semble etre le cas par exemple 
(quand on fait des simulations numeriques) du quart de plan {^z > 0, > 0). 
On sait en tons cas que Yq existe sur une partie de I'axe reel, a partir d'un reel 

>0 jusqu'au voisinage de I'infini. 

2.3 Existence d'une solution de (E'), Y{X,e), qui est voisine 
deFo 

Nous utiliserons un theoreme de point fixe pour montrer que certaines solutions 
de E' tendent vers Yq quand e tend vers 0. 

(E') YY' = 2XY (l - ^X^ + a + 1 - e'X, avec e' = £^/^ 

On ecrit Y sous la forme Y^Yo + e'Z. On rappelle (2.2) : YqY^ = 2XYo + 2. 
Alors 

YoY^ + e'Z'{Yo + e'Z) + e'Z {2X + |- 

= 2X^0 - s'X'^Yq + e'ZX{2 - e'X) +a-l + 2-e'X 

^ Z' = -X^ + 



.2 , a/e'-X-2Z/Yo 



Yo + e'Z 
avec a — a — 1 — 0{e) — o{e') . 



Apres linearisation de cette equation differentielle, on obtient 

Z' ^ - X^ + uo (^^ - - z(ul{2 + a - e'X)^ + e'h{Z, X, e') (2.5) 



et h{Z, X, e') = Z^ul 



On introduit les fonctions 



2 3 2 + a - e'X 



1 + e'uoZ 



(2.6) 



F{x) = J -ul{t){2 + a - e't)dt, 
et R{x) = ^F{x)) . 
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On designera par 7x un chemin dans le plan complexe partant d'un point 
(eventiiellement infini) a determiner et allant jusqu'a X. 

En integrant I'equation difFerentielle avec la methode de variation de la con- 
stante, on voit que Z satisfait a I'egalite : 



Z{X) = e^(^)-^W (^-t' + uo{t) {^-t)+ e'h{Z, t, e')^ 



dt. 



On souhaite montrer qu'on pent appliquer un theoreme de point fixe, pour 
montrer qu'il existe une et une seule solution a I'equation (E') qui tend vers Yq 
quand s tend vers 0. On se place pour cela dans un domaine ouvert V (qui sera 
precise par apres) et on considere deux solutions notees Yi et Y2 arbitrairement 
proches de Yq. 

Soit alors hi = h{Yi,t,e'), pour i — 1,2. 

On trouve 



Uo{t){2 + a-e't) uo{t){Yi-Y2) 



{Yi + Y2)uo + e'Y,Y2ul 



On est done amene a introduire les espaces fonctionnels suivants : 

Sj — Ifonctions holomorphes h sur V / sup 1 — , . , ^^^^^^ soit fini^ 

I- / tGv\uo{t)\-\2 + a-s't\ J 

2) = < fonctions holomorphes Y sur T> / sup \ Y{t) \ ■ |mo(^)| soit fini > 
1^ / tev J 

que Ton munit des normes correspondantes \\h\\sj et \\Y\\<q respectivement. 



Soit Oh : 2)^io et Oy : ^ ^2) . 

Y^h{Y,t,e') h^Z{X) 

Ces deux operateurs sont bien definis an moins sur des sous-ensembles Jq et 2) 
des espaces et 2) (ceci sera reprecise plus loin). 

Nous allons montrer que I'operateur OyoOh est une contraction sur ces ensem- 
bles. 



Il^ilb 


+ 


ll^2||a) + 


e' 


■ ||yi||2)||^2||2) 


(1- 


e' 


■\\Yih){l- 




e' 


• ll^2||sj)) 



\\hi-h2U<\\Y^-Y2\\^ 
done si ||Yi||a) et ||5^||2) sont majorees par 5, et pour e' < 1/5, on a 

\\h, - h^U < \\yi - y2h n LL2 ■ (2-7) 
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Par ailleurs, si Yk{X,e') — Oy(/ifc(X, s')), 
Y,-Y2^ f e^(^)-^W£'(/ii - h2)dt 



d'ou 



\Yi-Y2\\^< sup \uo{X)\\e'\ [ e^(^)-^(*)|/ii - /isl^i^ 



< \6'\.\\h-h2hsnp e^w/ e-^(*) (2 + a - e't) 



Mo(X) 



uo(i) 



Choisissons un X tel que : pour tout X' tel que \X'\ > \X\, uq{X') est peu 
different de —X'. Plus precisement, on va choisir 7x et X tels que : 

1- 7x(^) £ 1^ pour tout t. 

2. 7x va de I'infini a X avec jxit) — e'^i, ou t est un reel allant de +oo a \X\, 
et 9 — argX. 

Alors, pour 9 dans I'intervalle j ^ + f ~ <^ , 



F{X)-F{t)^ -ul{(l)){2 + a-e'(f))d(f) t ^ \t\e 

J\X\ei» 

r\t\ f\t\ 
= / (2 + a)(-0')d(/) + £' / (tf'd(t) + o{\t\^) 
J\x\ J\x\ 



La partie reelle de cette derniere expression est positive, (pour \t\ assez grand 
donne independamment de pour tout arg(£:') assez voisin de 0) et croissante 
avec \t\. Done si on choisit un X assez grand, il existe un // > tel que 



jR{lx{t))\>\F'{^x{t)) 



, pour tout t, 



(2.8) 



d'ou on deduit 

\\Yi - Y2\\<Q < \s'\.\\h, - h2h sup e«(^) / 

X€V Jj 

\\Yi - Y2\\<x) < \s'\ ■ n ■ \\hi - /islU 



dt 



(2.9) 
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Nous pouvons maintenant preciser le domaine T> : 

• V C (Domaine de Iq) (pour que Uq soit definie). 

• Tout point de X e "D pent etre joint de I'infini par un chemin 7x verifi- 
ant (2.8). 

• (2 + a — e't) ne s'annule pas pour t dans ce domaine. 

Le domaine T> n'est pas vide : avec les conditions donnees, on pent prendre 

par exemple pour T> un secteur b I oo, 2|A/|, + 77 ~ " J ■ 

V 8 2 y 

Alors, pour e assez petit, les inegalites (2.7) et (2.9) signifient que OyoOh pent 
etre un operateur contractant. 

II ne reste qu'a verifier qu'on a bien Oy(i5) C 2) et que Oij(t}) C ou t) est 
une boule ouverte de 2). 

Si y e 2) et ||V||2) < — -f (on a la la description de I'ensemble X)), alors 
Oh{Y) = h{Y, X, e') = uo{X){Y{X)uo{X) f 2 + a - e X 



\\h{Y,X,e')\\s^<snp 
xeT> 



1 + e'uo{X)Y{X) 



1 + s'uou 



qui est bien borne, et a Y fixe est aussi borne pour tout e' assez petit, puisque 
\uoY\ est borne pour X & V quand F e 2). 
Si heS^, 

OY{h) ^ Z{h,X,e') 

^F(x)-F(t) l^_^2 ^ - ^) + e'h{t, s')j dt 

= Z(0, X, e') + (Z(h, X, e') - Z(0, X, s')) 



La diff'erence dans la parenthese ne posant pas de problemes vu ce qui a ete montre 
ci-dessus (la fonction nulle est clairement dans I'espace i^), on regarde uniquement 
le premier terme. 



Z(0, X, £') = / e^(^)-^W (-t^ + u^it) {-,-t\) dt 
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On retroiive quelque chose d'analogue a la demonstration concernant la difference 
des fonctions Yi —Y2. 

\\Z{0,X,e')\\y< sup \uo{X)\ I e^(^)-^(*M-tV + Mo(t)(a-£'t)Ut 

< J_. sup e^(^/ e-«W \-euoe' + (a - e't)' 
|£ I Xev 

\\Z{0,X,e')\\^<\\-uoe't{uo{t)+t)+aul{t)\\^supe'^^''^ I 

xev Jj^ 



Uo{t) 



dt 



dt 



R{7x{t)) 



dt 



La fonction —UQs't{uo{t) +t) + aul{t) est bien dans Sj : uo{t) — —t + o{l/t), done 
la fonction est majoree par \uo{t){l + \at\ + \s'\)\ < \uo{t){2 + a — e't)\ pour tout 
e assez petit. 



\\Z{0,X,e')\\fQ existe bien, Z{0,X,e') e ^r) . 



On pent done appliquer le theoreme du point fixe, ce qui nous donne I'existence 
et I'unicite de la solution de (2.5), Z{X,s'), dans I'espace 2) : cette fonction est 
holomorphe et tend vers quand X — > 00, pour tout e'. 

On en deduit : 

Theoreme 4 II existe une solution de (E') sur un domaine V contenant un 
secteur ouvert centre en l'inji,ni (de taille independante de \e'\) d'ouverture un 
peu inferieure d 7r/8. tend vers en I'infini et elle est telle que Y^ — > Yq 

uniformernent. 

Ce domaine T> peut meme etre encore etendu. 



Montrons qu'on peut, en effet, prendre un domaine T> plus grand que celui 
propose. Comme I'equation (E') n'est pas singuliere, les chemins peuvent etre 
rallonges d'une distance arbitraire (mais finie) a partir de Xf (par exemple, on 
peut aller jusqu'a 0, ou Xi) ; I'integrale 



uo{X) 



Uo{t) 



dt 



restera effectivement bornee si on evite les singularites de Yq, puisque toutes les 
fonctions qu'elle contient sont holomorphes et done bornees sur tout compact 
inclus dans leurs domaines de definition. Et puisqu'il n'y a pas dans cette integrale 
de singularite en e, on pourra choisir des bornes independantes de e. 
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On pent menie etendre T> sur une partie de I'axe reel, et ce jusque ' 
par exemple. En effet, supposons que Ton ait deja prolonge 7x de X a Xi, puis 
eventuellement a G M"*", Xq assez grand. 

Alors 

R{x) = Cste + / -t^{2 + a - e't)dt + o{x^) 
Jxo 

^Rit) = - t^i2 + a - e't) + oit) e R" si t<- 
at e' 

(pour s' assez petit, a etant de I'ordre de e'^) 

\\Yi - Fall?) < \e'\.\\h - /i2|U sup e«(-) j^'^' e'^^ 

< - /i2|U Cste 

ce qui implique une convergence uniforme sur puisque I'inegalite (2.7) 

reste bien sur verifiee tant que Ton ne s'approche pas des singularites ou des zeros 
de la fonction uq. 

Cette solution de (E') correspond a une solution v de (E) qui existe done, 
pour £ > 0, au moins jusqu'en 2|X/|e*^'^/^^ x e^^^, sur le segment [—1 + Xis^^^, 0] 
et qui coincide au voisinage de I'infini avec la solution y etudiee dans [7] : ce sont 
des prolongements holomorphes I'une de I'autre. 

Cette solution admet un developpement asymptotique uniforme : 

Corollaire 5 // existe une suite de fonctions Yn{X) holomorphes en X et tendant 
vers quand X tend vers Vinfini telle que, pour tout N & N, et X &T>, 

N 

Y,{X,e')-Y,yn{x)e' 

n=0 

oil la fonction bornee verifie en outre 

lim Rn{X) =0. 

Le cas N — correspond au theoreme precedent. On va demontrer le cas 
general par recurrence, a partir du cas N — 1. On reprend I'equation 2.5 : 

Z' = -X^ + uo (^^ - X^ - z(ul{2 + a - s'X)^ + e'h{Z, X, e') . 

On a montre que Z(X,e') etait une fonction qui existait, dans I'espace 2), et que 
sa norme ||-Z^||a) etait bornee independamment de 



<£'^+^|i?;v(X)|, 
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On definit Yi comme solution de I'equation differentielle ci-dessus avec e' — : 

Y;{X) = -X' - XuoiX) - 2Y,{X)ul{X) , 

bornee dans V. Dans le secteur [3tt/8 + 6, 11/2 — 6] considere, le coefficient du terme 
homogene, —Uq{X), a une partie reelle qui tend vers +00. II existe done une et une 
solution Yi de I'equation differentielle qui soit bornee dans le secteur considere. 
De plus, le terme non homogene est d'ordre au plus —X'^ — Xuo{X) — o(l) ; on 
peut en deduire, comme pour Yq, que Yi{X) tend vers en I'infini. 

On pose maintenant Z{X,e') = Yi{X) + e'Zi{X,e'). 

Alors Zi{X,e') est solution de 

£'Z[ = uo + s'hiZ, X, e') - (Fi + e'Zi) (ul{2 + a - e'X)^ + 2Yiul 
Z[ = -Z,{2 + a- e'X)ul + (^) + ulYi (^X - ^) + h{Z, X, e') , 

ce qui est une equation differentielle lineaire en Zi (le fait que Z puisse s'ecrire en 
fonction de Zi n'est pas genant : on salt que cette fonction Z existe, et h{Z, X, e') 
peut etre considere comme une donnee, comme fonction de S3). 
On peut resoudre cette equation : 

Z^iX, e') = e^(^'^') e-^(*'^') (uo{t) (^) + ul{t)Y,{t) (t - ^) + h{Z, t, e')) dt , 

en rappelant que F{X, e') — J —ul(t)(2 + a — £'t)dt. 

, / ,M h(Z{t,e'),t,e') 
Alois, puisque |'Uo(^)(2+a— ^ ^)l > dans V, et que les fonctions — 

Uo{t) 

et Yi(t)uo(t) sont bornees, il est clair que Zi{X,e') est une fonction bornee, qui 
reste dans I'espace 2). Ce qui demontre le cas = 1. 

Pour la suite de la demonstration par recurrence, elle devient evidente a par- 
tir de I'equation differentielle obtenue pour Zi : on va trouver a chaque fois 

- F„_i(X) + puis 

Z; = -Zi(2 + a - e'X)ul + e') , 

avec hn{X,e') — ^"■"^^^'^^^ ^hn-i{X,0) ^ 2ul{X)Yn-i{X) . hn est une fonction 

de 9) avec I'hypothese de recurrence qu'on peut exprimer comme ceci : /i„ G et 
Zn-i € 2). Done, en resolvant I'equation, on voit a nouveau que Zn{X,e') est une 
fonction de 2). 

Ce qui suffit a demontrer le theoreme. 
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2.4 Estimation de — a 
2.4.1 Preliminaires 

Premiers termes du developpement de v{u,e) et Y{X,e) 

Le calcul des premiers termes des series pour a et v (cf (2.1)) nous donne 

a{e) = aiE + 02^^ + 0{e^) v{u, e) = vo{u) + evi{u) + e^V2{u) + 0{e^) 



ou 



1 3 

tti = — - 02 = — 



et 
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, , -1 , , 1^2 + 4m + 7 



1 3u^ + 2\u^ + 66?/ + 126^2 + 159k + 121 

= -32 (^TTIF 

Par recurrence, on demontre que pour n > les f„(M) sont des fractions ra- 
tionnelles, dont le numerateur est un polynome de degre 3n — 1, et dont le denom- 
inateur est {u + 1)^"+^ 

Le simple changement de variable u——\ + e^l^X dans le developpement de 
v{u^e) nous donne alors, apres simplification : 



X 2X^ AX^ \AX^ 8X^ J \8X^ AX^ 



SX^ 16X3 32X2 \^|J5C|8 I l|J^|5 I I \x 

Ensuite, comme Y{X,s) = s^^^v se met sous la forme 

Y^Yo + e^/^Yi + ■■■ + y„£"/^ + o , 



35 



apres identification, on trouve pour les : 





1 

X 


Yi^ 


1 


4X3 


Y2^ 


1 


8X2 






et 


pour n - 







2X4 ^ \^JSC7 

V-^V (2.10) 



Equivalent de Yq' — Yq 

On introduit les fonctions suivantes : 

Aio = Ai, Ail : X (-^ Ai(ja;), Ai2 : x i— > Ai(j^a;) . 

On cherche a estimer la difference de deux des solutions de (2.2). Pour cela, nous 
nous interessons aux deux fonctions suivantes : 



Aii(.-) " ' Ai2(.-) 

Ces deux fonctions verifient, a ca use de la symetrie par rapport a I'axe reel de 
I'equation (2.3), la relation Xi{z) — X2(^). Elles correspondent a deux solutions 
de I'equation (2.2) : 

yo+(X) = Xl {z{X)) + z{X) et Y,-{X) = Xl {z{X)) + z{X) . 

En particulier, pour Xi[z) sur I'axe reel, on a : 

X :=Xi(^) = X^{y,- {X,{z)) - Xl{z)) = X,{y^-{X) - X^) 

= X, {z) = X, {X) - X2) = X, ( V {X) - X^) . 

On note dans ce cas 

z+{X) =Y+{X) - X^ et z-{X) = Fo"(X) - X^ 

On cherche un equivalent de Y^^ — Fq", et nous allons montrer le resultat 
suivant : 
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Lemme 6 Pour X sur I'axe reel, 



x^+oo e 

Ce qui revient a reussir a calculer exactement un coefficient de Stokes pour une 
equation difi^erentielle non lineaire ; dans le cas general, ce n'est pas possible. 

Pour obtenir ce resultat, on utilise le fait que pour X reel et suffisamment 
grand, FJ^ — X'^ est presque un reel negatif. Dans ce cas |Aifc(y^'^ — X"^)] est grand 
pour k = 1 eX k = 2, exponentiellement petit par contre pour A; = 0, ce qui 
justifie I'utilisation d'un developpement limite pour estimer les deux fonctions 
A\2{z) = —j^A\i{z) — jAio(^) et Ai2(^) = —j^A\[{z) — jA\Q{z). Dans le calcul qui 
suit, on travaillera en variable X, avec aussi les fonctions z^{X) introduites plus 
haut. 



X, (Fo-(X) - X') = X2 {Y,+ {X) - X') 



Ah{z-) Ah{z+) 
p AiU^-) ^ Ai[{z^)+fA\',(z^) 
Ah{z-) Ah{z+)+fA\o{z+) 

II va etre possible de calculer un equivalent de cette difference, ce qui nous menera 
a un equivalent pour — Fq". En effet, avec une formule de Taylor 

X,{z-) - X,{z^) r^iz- - z+)X[{z-) = - Y,+)X[{z-) 

sachant que d'apres I'equation differentielle pour Xi, on pent ecrire sa derivee 

Xf{z-)+z- 



X[{z-) 



X[[z-iX))^ 



2 

+ 1 



2X 



Nous allons utiliser les equivalents deja vus precedemment (cf (2.4)) pour 

Ai' (z+) Aio(-2+) 

estimer les deux fractions — — et — — — — , puisque z^ est proche du demi-axe 

A\^{z+) A\i{z+) 
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reel negatif quand X — > oo : 

=i(x»4+o(i)) 



3/2 



Ai'o(z+) Aio(z+) _ -e-«zV4 ^ 1 _ ^ X ^^^^^^^ 



En rappelant que 

on trouve comme equivalent pour la difference le resultat annonce au lemme 6 : 

(Y+ -Yq-)(X) r^-X^e-i^\ quand X ^ oo dans K. (2.11) 
e 

2.4.2 Etablissement de I'equation verifiee par a'^ — a" 

On rappelle I'equation (E), verifiee pour {v,a) — {v'^,a'^) et {v~,a''). 

dv 9, a — u 

E u') + . 

au V 

Le but des paragraphes a venir va etre d'obtenir une estimation de — a~ . 
Pour cela, on commence par chercher une equation differentielle en — v~, ou 
le terme a calculer sera present. On pose alors 

b=a~^ — a~ c=a^ + a~ 
d = — v~ e—v^ + v~ 
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Quand on effectue la substitution dans I'equation (E), on obtient 

ed' = (1 - u^) + - (1 - u^) - 

ud + a'^v~ — a~v'^ 

b+c e—d I b—c e+d 



M , e , c , 
= — — d H — h - — - — d 

v~^v~ 2v~^v~ 2v'^v~ 

u — c/2 , + v~ , 

= , b 

v~^v~ 2v~^v~ 

En remplagant c = + a~ = 1 + a"*" + 1 + a~, on obtient en definitive 
I'equation voulue. 

Lemme 7 On pose d{u) = v^{u) — v^{u) et b = — a~ . Alors V equation 
differentielle suivante est verifiee : 

ed'{u) = f{u)d{u) + g{u)h 



avec 



sachant que 

a = ^--e + 0{e^). 



(toutes les fonctions ici dependent aussi dee, cette dependance n' est pas rappelee). 

On integre ensuite cette equation diff'erentielle avec la methode de la variation 
de la constante, entre +oo et x (car lim+oo d = 0) : 

^ J+oo 

ou 

F{t)- fin)du. 
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Le choix de x — xi — —1 + XiS^^^ permet de trouver que 

h = £d(x,)e-^(-')/- . (2.12) 

2.4.3 Calcul de I'equivalent pour b 

On peut trouver une estimation, quand £ — > 0, de chacun des facteurs dans le 
produit de I'egalite (2.12) ; le but final etant d'en obtenir une pour b. 

Tout d'abord, on a vu que jusqu'au point u = Xi, on pouvait aflirmer que la 
fonction v+{u) = e~'^/^Y+{s~'^/^{l + u)) etait equivalente a e~'^/^Y^{e~'^/^{l + u)) 
quand £ — > 0. Ce qui nous mene a 

d{xi) = v+{xi) - v-{xi) ~ £-V3 fY+{Xi) - Yo-{Xi)) . (2.13) 

£—♦0 \ / 



Ensuite, I'integrale presente dans (2.12) est susceptible d'etre estimee avec la 
methode du point col. On constate d'abord que 

F'{t) = f{t) = 0^ t=l + a = a. 

Or on verifie sans difficulte que F{a) — 0{e'^). D'oii 

£ 2£ £ 

Sachant que i;^(q;) — —\ + C'(£), on trouve 



v~^{a)v (a) {v~^{a)v (a))" 

1 



-1/2)2 

4 + 0(£), 



+ 0{e) +0 



done 



Par ailleurs 



F{t) _ 2{t - ^ 



g{a) 2. 



2 0( + k|.k-a'^ 



) 
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La contribution principale dans I'integrale vient du voisinage de a, on va done 
ecrire, apres le changement de variable x = {t — a)^/2/e, 




Un equivalent de I'integrale est done : 

r+oo 

- (7(t)e-^WAdt~y2^. (2.14) 

J XI 



II reste a trouver une estimation pour F[xi). Pour calculer F(xi), on pent 
couper I'intervalle [xi, 1] en deux parties, en un point (—1 + A), ou A est un reel 
fixe a choisir entre et 2. A gauche, on pourra alors raisonnablement estimer v{u) 
par I'expression e-^/^'{Yo{X)+e^/^Yi{X)+e^/^Y2{X) + 0{e)) d'apres le corollaire 
5 , et a droite par Vq + evi + 0{e'^). 



1 , 1 /■"'^^ -u-l-a , 1 r 
-F{xi) = - / — -- — --du+- 
e e Ji v~^[u)v^[u) e J_ 



(A) (B) 

Les calculs des premieres fonctions Vn et les approximations pour les F„ (cf. (2.10)) 
donnees dans les preliminaires permettent de donner une estimation des deux in- 
tegrales. Pour I'expression (A), on part de 



(^ + 1)1-7 , .V. +^(^ 



v^{u)v {u) \ 4 (it + 1)^ 

puisque v"*" et v~ ont le meme developpement asymptotique, d'ou 



e Ji v+{u)v [u) 
= -/ (u-l-a)(u + lYdu- - (u-l-a) du + 0(e). 
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On rappelle que a — —e/8 + 0{e'^). Le calcul des integrales nous donne alors 



Pour (5), le calcul se revele plus complexe. On salt que chacune des deux 
fonctions v'^ et v~ se prolonge a I'aide d'une fonction du type du theoreme 4 : 

avec X — u/e' , en utilisant le corollaire 5, 

ou le Os{l) est borne quand s ^ et au moins borne quand X 00. 

Comme et v~ out le meme developpement asymptotique pour tout u > p 
dans un certain secteur autour de I'axe reel, il est clair que le developpement 
en 1/X des Yn donne plus haut est valable pour Y^ comme pour Y~. Ces deux 
fonctions sont done exponentiellement proches en 1/X quand X devient grand. 
On pent done ecrire, 

IT' M - 1 - a , 
S = - / ^-——-du 

e J_i+y,v+{u)v iu) 

en notant 

_ Y+ yf _ y+yf yf ^ y2+ yg" 

le reste R{X, e) se trouve etre de la forme + e^/^^ ^e{^), puisqu'il correspond 

y± ^4/3 ;^ 

a des termes (1 < /c < 3) et —7^77^ respectivement equivalents quand X 



'0^ -"O 

devient grand a — et e^^^. 

X 



A.-1/3 yo+yo" ^ A.-V3 yo+yQ-^ 
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On appellera les quatre integrales precedentes, dans I'ordre, (/), {II), (HI) et 
{IV). 



Pour (/), regardons d'abord f, L1/3 +?,_ dX. D'apres I'equation differentielle 
(2.2) verifiee par Yq, on a 

-2 

Y,+{X)Y^{Xy ^ 

=^ (V - = (V - ^o")(0)e (en integrant I'equation) 

Or, quand X — > 00, on a vu avec (2.11) que 

4 .„2„-|X3 



{Y+-YQ-){X)r.-iX'e-s 



done 

A.-./, ii+ro' V i / V ' J 

= ln(-i(yo+ - Yq-)(Xi)) + 1 - 2 ln2 + — - 2 ln(A£-V^) + o(l) 

= ln(-z(Fo+ - Y,-){Xi)) + _- + -Ine - 21n2A + 1 + o(l) . 

Pour tons les termes a venir (y compris pour les autres integrales), I'integration 
au voisinage de Xi et de tout point fini ne donne au final que du o{e), et il n'y 
a done que ce qui se passe quand X ^ 00 qui nous interesse. Autrement dit, 
la difference entre Y^ et son developpement en 1/X est negligee. Ainsi, pour 
fx -1/3 v^i^dX, on utilise le developpement asymptotique de Yq : 



pXi 1/3 Y r 

/ ^dX^s^'^ {X'-l)dX + 0{e 



l/3> 



£ 4 



Et pour la derniere partie de (/), 

-dx(l + 0{e)^ 



X3 
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D'ou finalement 



1 /2A^ 2 

(/) = -(— -) + -ln£-21nA-21n2 + l + A 

£ \ 3 4 / 3 



On peut proceder de maniere tout a fait identique pour (77) 

en utilisant la serie asymtotique de Yi dans Qi et ne negligeant ce qui est en 
facteur de a = 



{II) = e 



1 



+ o{x- 



1/3 



4X3 

^="'^'-2 + £V3X 



21^ 



-2+0(1) 



A^ 



dX + o{l) = — - A + o(l) 



Et de meme pour (III)- 

(///) = "(-2 + e'/3x) (_^:) = ^ _ ^ + 



Pour la derniere integrale, 

Ja£-1/3 

<£C»(£-^/^) 
= 0(1) 
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En fin de compte, quand on additionne (A) + (I) + (77) + (777) + (71^), on 
obtient 



£ £3 3 3 \ I / ' ' 



Done 



-F(xi)/e ^ g-4/3£ -2/3^^ I (2 1 5^ 

Remarque 

On peut signalor quo lo calcul dc (A) ci-dessus n'cst pas indispensable et qu'il sufRt dc poser 
A = dans (B) pour obtenir le bon resultat. Le calcul complet de {A) a cependant ete conserve 
pares qu'a priori, il est indispensable, qu'il peut I'gtre dans d'autres exemples et parce qu'il 
permet de detecter d'eventuelles erreurs. Dans I'autre sens, remplacer A par Xie^^^ dans {A) ne 
sufRt pas pour obtenir le resultat complet, mais donne deja les termes e~'*/'^^e~^/^. 

On peut en revanche choisir effectivement n'importe quel A fixe dans ]0, 2[, ou m6me un A 
variant tres peu avec e comme par exemple e^^^^, on aboutira au bon resultat, avec quelques 
complications dans le second cas. 

En combinant les resultats intermediaires (2.13), (2.14) et (2.15), eela nous 
donne finalement : 

^ ^ e4/3 (F+-r-)(X,) 

b ~z 



Soit le 

Theoreme 8 



a —a ~ z 



On retrouve le resultat de [7] de la proximite exponentielle des deux a. Et on 
confirme que la borne trouvee : exp (7?(— 1 + p)/e) etait bien optimale. 

2.5 Consequence sur les coefficients de la serie k{e) 
D'apres [7], la fonction a{e) est holomorphe dans le domaine 



X = <( £/arg(£) e ^ + 5, 1 -(5 



,0 < |£| < |£i| >. 
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Pour £o assez petit (|£o| < le chemin 7, defini par un arc de cercle de rayon 
|eo| et deux segments proches de [O.^o] (voir figure 6) est inclu dans A4, et est 
retractile en un point dans ce domaine.. 




Fig. 6: Chemin d'integration 7. 

Done pour e G on a 

a(u) 



a(e) = / 



u — e 



-du 



1 t -M^^^^r(a(ne-')-a{a) 



2m J\u\=p u — s 2m Jo V u — e 
puis en developpant en serie de e et en identifiant les coefficients obtenus 

a^ = — f a{u)u-''-^du + — f ' ( a{ue-^''') - a{u)]u-''-^du 
^^T^ J\u\=p 2m Jq \ J 

comme le premier terme n'est que de I'ordre de p"" 

a„ ~ / {a~{u) — a^{u)\u'^~^du 

2m Jq \ / 

et que dans cette integrale, la contribution principale vient du voisinage de 0, le 
theoreme 8 mene a 

1 . 8\/2 f'0^-A/3u 
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On effectue alors le changement de variable t — -^,du — — |^ 

Pour n assez grand, /"^ e^H^'-^l'^dt est negligeable face a e-H^'-^/'^dt. 



Done 




Theoreme 9 On a un equivalent exact pour les coefficients des and canards » ; 

7re4/3 {AeJ 

On rappelle que tons les a donnant des solutions « canards » sont exponentielle- 
ment proches en e. lis ont done tons le meme developpement asymptotique. 
Par ailleurs, le resultat obtenu (a e donne) en sommant la serie «au plus pe- 
tit terme » donne aussi une approximation de a exponentiellement bonne, et est 
done une valeur a canard [7]. 



2.6 Resultats numeriques 

On pent facilement implementer le calcul de la suite des a„ en utilisant les 
formules de recurrence (2.1). En particulier, en calculant les 155 premiers termes 
(calcul effectue par Franck Michel et amicalement transmis) et en regardant le 
resultat de la multiplication 6„ = a„ x (4e/3n)", on obtient : 
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&135 


-0,5417512651 


^136 


-0,5418690317 


^137 


-0,5419854885 


^138 


-0, 5421006603 


^139 


-0,5422145711 


^140 


-0, 5423272443 


bui 


-0, 5424387024 


bl42 


-0, 5425489682 


bl43 


-0,5426580621 


6l44 


-0, 5427660064 


bl45 


-0, 5428728208 


bl46 


-0, 5429785257 


&147 


-0, 5430831405 


bl48 


-0, 5431866841 


bl49 


-0, 5432891757 


^150 


-0, 5433906324 


bi5i 


-0, 5434910728 


bl52 


-0,5435905137 


&153 


-0, 5436889722 


bl54 


-0, 5437864645 


bl55 


-0, 5438830066 



Si on cherche alors les 6„ sous la forme bn — C + ajs/n et bn — C-\-a/ \/n, on 
obtient (en utilisant une formule de linearisation aux moindres carres) : 



bn = -0, 5736898877 + 



0, 3710889332 



bn = -0, 5891153498 + 



0, 2429519729 



n 



Et la constante calculee ci-dessus vaut 



-0, 5813148764. 

On constate done que la constante calculee est compatible avec les resultats 
numeriques, mais que les fitermes correctifsz du developpement des a„ restent 
non negligeables. 



II est aussi possible d'estimer directement la valeur deb — — a~ — 2'is{a'^) 
en cherchant les valeurs numeriques de a'^ pour differentes valeurs de e. 

On integre numeriquement, pour certains e raisonnablement petits, I'equation 
differentielle (E), par example d'une part le long du chemin comprenant les seg- 
ments [— 1 + lOi, 0] puis [0, 1] et d'autre part le long de [9, 1]. Ces deux chemins ont 
la propriete de descendre le relief (cf. fig. 7, et [7]), c'est a dire qu'une erreur dans 
I'estimation des valeurs initiales (en — 1 + lOi et en 9) devient exponentionnelle- 
ment petite une fois arrives en 1 ; on prendra done les valeurs f ( — 1 + lOi) = i/10 
et v{%) — —1/10. De plus, un de ces chemins part de la montagne Nord, I'autre 
de la montagne Est, et la definition de est qu'il s'agit de la valeur de a pour 
laquelle il existe une solution tendant vers en Too sur ces deux montagnes (N.B. : 
on retrouvera la maniere de construire les a dans la partie suivante), ce qui sera 
vrai si et seulement si les deux v{l) trouves sont egaux. 
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Fig. 7: Relief pour I'equation de Van der Pol, et chemins pour I'integration numerique 
Une recherche numerique de ce type donne les approximations suivantes : 



e 


«+ 


2$j(a+) X e^/3^V^ 


0,20 


0, 9684 + 0, 00153i 


1,07 


0. 17 


0. 9733 + 0.00055/ 


1.16 


0,14 


0,9800 + 0.000120^ 


1,23 


0,08 


0,9893 + 1,40- 10-''^ 


1,37 


0,06 


0,9921 + 6,48- 10-^^i 


1,42 


0,05 


0,9935 + 8,5- 10-^^i 


1,44 


0.01 


0.9918 + 1.23 • 10""/ 


1.47 



Ces resultats sont compatibles avec la valeur de la constante trouvee (egale 
environ a 1,68), avec des termes correctifs en ^/e et \/e. 
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Chapitre 3 

Generalisation 



3 . 1 Preliminaires 

3.1.1 Rappel : estimation de Cauchy 

On rappelle pour les demonstrations a venir la propriete de Cauchy, qu'on 
utilisera sous la forme suivante : 

Lemme 1 Soitf{xi, . . . ,Xn) une fonction holomorphe enx = {xi, . . . ,Xn), bornee 
pour X dans Vouvert U — Ui x U2 x . . . x Un), et soit Fi un ferme indus dans 
I'ouvert Ui. Alors 

df 

— — (x) est bornee pour x & Fi x U2 x . . . x Un 
0x1 

Un corollaire immediat de ce lemme indique que si ce n'est pas f{x), mais 
4^ qui est bornee, alors la fonction -i^S^(x) est bornee elle aussi. 

3.1.2 Forme normale pour une equation differentielle 

Plutot que travailler sur une equation ecrite sous la forme (E.diff.) vue dans 
I'introduction : 

^y' = yfi^, £) + h{^, £) + ^y^Pix, e, y) 

nous allons manipuler dans cette section une equation sous une forme (l) ap- 
paremment moins generale : 

ey' = {x"f{x) + eg{x, e))y + h{x, e) + ey^P{x, e, ey) . (l) 

En fait, on se ramene facilement a une equation de ce type. 
Soit I'equation 

eu' = ^(x, u{x),e) , 
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ou on demande juste que ^ soit holomorphe en x, meromorphe en u et qu'elle 
admette un developpement asymptotique en e quand e tend vers 0. On suppose 
aussi qu'il existe une courbe lente uq{x), telle que ^'(a;, Mo(a;), 0) = 0. 

On pose u{x) = Uq{x) + ey{x), ou, ce qui donne souvent un resultat plus 
agreable a manipuler, u{x) — Uo{x)(l + ey{x)), et on linearise I'equation en y. 

eu' — ^{x, u, e) , 

donne, en regroupant dans ipo et ipi tout ce qui est lineaire en u, et dans le dernier 
terme tout ce qui contient un pole en u — 0, 

= ipo{x,e) + uipi{x,e) + u^ip2{x,e,u) + ■\ip3{x,e,u) 
On ecrit alors u = Uo{l + sy) : 
£1^0 (l + ey{x)) + e^uoy' = V'o + ^^oV'i + ^uoy 

+ uli^2 + ul{2ey + e^y'^)ip2{x, e, uq{\ + ey)) + ^ \ — ^^3(2:, £, iio(l + ey)) 

Uq\i + zy) 

Comme uq est courbe lente, les termes de degre en £ s'eliminent ; on divise alors 
par £, puis on linearise autour de y = : les termes en y*^. A; > ont done comme 

z^ 

facteur un terme — , ce qui donne bien en particulier, pour /c > 2, £y'^P(x,s,£y). 

On obtient done bien une equation du type annonce ci-dessus, sachant que les 
fonctions f,g,hetP contiennent eventuellement des termes avec des facteurs en 



Ainsi, dans I'exemple de I'equation de Van der Pol (E), cette procedure de 
normalisation nous donne : 

(1 — u'^)v + a — u 

(1 - ^.^) + . 

V 



ev 



dv 

du 
dv 

du 



51 



Puis en posant v — vq{1-\- ey), vq etant egal a 



a — u 



e{vo{l + ey))' = l-u' + 

VQ{l + ey) 

/ /N , 2 l-\-a-uf^ e^y"^ \ . . 
evo{ey ) H [l-ey + — - evo{l + ey) , 

ey' = y (-{1 -u){l + u)^ + e (-{1 + u)^ ^ 



e 1 + u 

+ -(1 + u)^ - + ey^-^ ((1 -u){l + u)^ + a(l + u)^) 

e^ ^ 1 + u 1 + ey ' ^ ' ' 

Pour cette equation, I'etude a ete faite au voisinage des deux points cols du relief 
construit a partir de — (1 — + uf' : autour de u = \, dont on a montre qu'il 
n'etait pas un point tournant dans certains cas ([7]) ; et au voisinage de —1, au 
cliapitre 2. Ce sont ces deux etudes que Ton souhaite generaliser ici. 

3.2 Existence de solutions d'une equation difFerentielle pres 
d'un point tournant 

Dans les paragraphes 3.2.1 et 3.2.2, nous allons generaliser le theoreme 4 qui 
a ete demontre dans le cas de I'equation de Van der Pol, au cas suivant : 

ey' = {x^j{x) + eg{x, e))y + h{x, e) + ey^P{x, e, ey) . (i) 

II est clair que est un point col pour le relief 3f? {^f'^ f{t)dt) correspondant, et 
c'est done en general un point tournant pour I'equation differentielle. On veut 
etudier le comportement des solutions de cette equation au voisinage de ce point. 
Le but va etre de demontrer que si on a les bonnes hypotheses sur les fonctions 
intervenant dans I'equation, on a I'existence d'une vraie solution bornee, d'abord 
pour des x tels que > p, ou p > est une constante quelconque fixee in- 
dependante de e (resultat classique, redemontre au paragraphe 3.2.1 ci-dessous, 
pour e assez petit). Ensuite, que cette solution pent bien souvent etre prolongee 
dans des domaines de la meme forme, mais avec un p tendant cette fois vers 
avec £ : |a;| > p{e) = Xie^ (Xi constante assez grande, mais independante de e, 
et r un rationnel strictement positif donne). Ce deuxieme resultat sera demontre 
dans un second temps, dans le paragraphe 3.2.2, et on verra que si la solution 
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prolongee n'est pas necessairement bornee (quand £ — > 0) dans tout son domaine 
d'existence, sa croissance reste controlee et n'est pas, en tons cas, exponentielle. 

Chacun de ces deux resultats ne sera en general vrai que dans certains do- 
maines determines non seulement par la norme de x, mais aussi par le relief; de 
fait, ces domaines correspondront a pen pres a des secteurs centres en 0, avec des 
X verifiant la condition donnee ci-dessus : \x\ > p{s). 

3.2.1 Existence d'une solution nloinz de 

On se fixe done une constante p positive, arbitrairement petite. 

On suppose qu'on a, pour I'equation (l), les hypotheses (H«) suivantes : 

• / est holomorphe dans Dq, ou Dq est un domaine ouvert, fini on infini, 
contenant 0, et / ne s'annule pas sur Dq. 

• g, h et P sont holomorphes en x sur Dq, sauf eventuellement en ou g ou 
h peuvent n'etre que meromorphes et avoir un pole. 

• g, h et P sont holomorphes en s dans des secteurs ouverts Sq, centres en 0, et 
bornees pour tout x G Dq fixe quand e — : on suppose qu'il existe q tel que 
les fonctions P, x^g{x,e) et x'^h{x,e) out un developpement asymptotique 
en s quand e tend vers 0. 

• P est holomorphe en ey dans un voisinage de 0. 

• On pose F{t) = j\uPf{u))du et R{t) = K (^^^). On considerera a ^ = 

arg£ fixe un domaine ferme D51 C Dq, tel que soit accessible a partir 
d'un sommet sq, qui est soit infini, soit dans Dq \ avec le relief R{t), et 
tel que Vx e D, |a;| > p; on reprendra, pour ce qui concerne I'accessibilite, 
les notations de la definition 1. Les domaines peuvent etre non bornes, 
si Do I'est. 

L'equation (l) a une solution formelle y{x,e) — Yl,'^=oyn{x)s"' pour \x\ > p, 
dont tons les coefficients yn{x) sont holomorphes en x. Cela se verifie en rem- 
plagant dans l'equation y par cette serie et en identifiant les coefficients de : on 
trouve alors des relations de recurrence donnant ?/.„ en fonction des yp et de leurs 
derivees (pour p < n). La premiere relation, obtenue en posant e = donne 

h{x,0) 

*<^' = 

Nous allons demontrer le theoreme suivant : 

Theoreme 10 On suppose que pour I'equation (l), les hypotheses (H») sont veri- 
fiees et que les fonctions 

1 1 9ix,e) h(x,e) 

xPf{x) xPf{x) \y\<s xPf{x) xPf{x) 



53 



sont bornees^ pour un certain S, x G Dq, > r (r constante arbitrairement petite) 
et s E Sq, \e\ < £o- On fait varier e dans un secteur Sq d'ouverture arbitrairement 
petite. 

Alors, pour tout reel positif p > r, il existe un domaine ferme eventuellement 
non borne T) C Dq, de la forme 

( F(x) \e\ ~| 
x\ > p, X accessible d'un sommet avec le relief^ > 

e&So ^ 

et une solution y{x,e) de (l) existant pour x G D qui soit bornee independamment 
de \e\ sur ce domaine. Cette fonction y{x,e) tend vers la fonction yo{x) quand £ 
tend vers dans Sq, uniformement pour x E T). 

Quelques remarques sur les hypotheses de ce theoreme : 

1. L'hypothese l/(^x^f{x)^ bornee implique bien entendu que / ne s'annule 
pas sur T>. 

2. Comme il existe un Xp tel que Xp\uPf{u)\ > \g{u,£)\,\/u E D et e assez petit, 
cela signifie que le (vrai) relief correspondant a 3? ^ i^^f{u) + £g{u., e)Jdiij 
ne serait qu'une petite perturbation du relief effectivement considere. 

3. Comme h{x, e)/ (^xPf{x)j est bornee pour tout e assez petit et element d'un 

secteur Sq, la fonction |/o(a;) est elle aussi bornee sur D. Ensuite, d'apres 
I'inegalite de Cauchy, y'o{x) fait aussi partie des fonctions bornees sur D. 

4. Les points et les reliefs ne dependent que de I'argument de e, et pas de 
sa norme. L'intersection se fait done uniquement sur le (petit) intervalle des 
arguments de e possibles pour 5*0. Comme le relief et les sommets dependent 
continument de arg£, en prenant 5*0 suffisamment peu ouvert, on est assure 
que le domaine T) sera non vide. 

Le theoreme revient a montrer qu'il existe une vraie solution holomorphe cor- 
respondant a la solution formelle donnee ci-dessus dans certains domaines D pour 
tout £ assez petit. On commencera par considerer le cas oil arg£ est une constante, 
avec le relief et le sommet s correspondants. On omettra le plus souvent dans la 
suite de rappeler la dependance de la fonction y en e. 

On cherche cette solution holomorphe sous la forme y — yo + L'equation 
(l) devient alors 

sz'{x) = {xPf{x) + sg{x, s))z{x) + yo{x)g{x, e) - y'^ix) + ^ h{x,0) 

+ {yo{x) + ez{x)fP{x, s, syo{x) + e^z{x)) 



^autrement dit, la partie lineaire homogene de l'equation difffirentielle ne devient jamais nfigligeable 
par rapport aux autres termes ou aux constantes 
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En utilisant la formule de variation de la constante, on obtient 

z{x) = \ I e^'^(y,{t)g{t,e)-y',{t) + ^^''^ ^ ''^^'^^K H{t,e,z{t)))dt 

avec 7a; un chemin partant du sommet s et descendant le relief jusqu'a x, et avec 
la fonction H{t, e, z) — szg{t, e) + {yo{t) + ezyP{t, e, eyo{t) + e^z). 

On souhaite appliquer un theoreme de point fixe, pour e assez petit fixe. On 
considere pour cela les deux operateurs suivants : 

7i : 3 C 3 ^io 

z 1-^ ^{z), 

la fonction 'H{z) etant telle que T-C{z){x,e) — H{x,e, z{x,e)) pour tout {x,e) e 
et 

H^z = U (y^ixUt, e) - .^(t) + ^^'^ ' ^^'^ + HH)) dt 

oil ^ et 3 sont les espaces fonctionnels suivants : 

3 = {z{x,6)/ z holomorphe et bornee dans T> x So} 

S) = I^H{x,s) I H est holomorphe dans D x Sq, 

H(x,e) , , . , ,1 

et sup — est borne mdependamment de £ > 

Ces deux espaces sont munis des normes correspondantes || ■ || et || ■ ||^. Avec 
les hypotheses du theoreme, les fonctions g{x,e), P{x,e,£5) et h{x,e) sont dans 
dP 

i^, ainsi que — — -, et avec une norme dans cet espace que Ton pent majorer 

d{ey) 

independamment de e. 

Alors 7i et Z sont des operateurs contractants. En effet, on pent ecrire (en 
omettant de rappeler les dependances de g, yo, z^ etHk en x et e), 



n{zi) - n{z2) = 



eg + (2£yo + e^{zi + Z2))P{x, £, eyo + e^Z2) 

Zl - Z2) 



,/2,o ^ 2 2^P{x^^^^yo + ^'^zi)- P{x,e,eyo + e^Z2 



{zi - Z2) 
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et on constate que, si \\zk\\ < S (on choisit ici 6 pour que ce soit vrai) et pour e 
assez petit, on se trouve devant une combinaison lineaire de fonctions de Sj, avec 
dP 

g, P et —, — -, multipliees par des fonctions bornees sur !D, telles que vq et 
d{ey) _ _ 

Z2. Cette combinaison lineaire est done bornee, et comme on pent mettre un e en 
facteur, on a : 



\\n{z^)-n{z2)\\^<\e\C\\z^-Z2\ 



(3.1) 



Pour I'autre operateur, 

\\Z(H^)-Z(H2)\\ < 



e e 



e'''^\\H,-H2\\^\t^m\dt 



en utilisant la condition d'accessibilite du domaine D, 

\Hi — H2\\^ r Rix)-R(t) 1 —d 



< 



e e R(t)dt 

1 r R(x)-R(t) 1 a; 
e e 



Done I'operateur Z est tel que 

\\Z{H,)-ZiH2)\\<M.\\H^-H2\\^ 
oil M pent bien etre choisi independant de e. 



(3.2) 



Verifions que les ensembles d'arrivee des deux operateurs sont bien ceux qui 
sont donnes. 
Si z(x,e) e 3, 



\\n{z){x,e)\\^ < \e\ ■ \\z\\ ■ \\g\\^ + ||i/o + ■ \\Pix, e, e6)\\^ , 



(3.3) 



ceci, si pour tous les z tels que \\yo + ez\\ < 6 (d'ailleurs, I'operateur n'est defini 
que pour de tels elements de 3, et non sur tout I'espace). 
Si H{x,e) eSj, 

Z{H) = Z(0) + {Z{H) - Zm 
\\Z{H)\\<\\Zm + \\Z{H)-Z 



1 

< - 

£ 



\Z{H)\\ < M 



yog -yo + -Q^ 



yog -yo + 



dh 
dh 



R{x) 



R(t) 



e-—\tPf{t)\dt + M\\H-0\\^ 



+ M\\H\\^ 



(3.4) 
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En effet, yog, z/'q et — sont dans I'espace S), puisque yo et y'g sont bornees, et que 
g e et h e Sj. 

Done, vu les inegalites (3.1) et (3.2) le theoreme du point fixe est applicable 
dans ce cas pour I'operateur ZoH pour tout e dans Sq, assez petit. 

Nous venons de demontrer exactement qu'il existe une solution y de I'equation 
dilferentielle etudiee qui tend uniformement (pour a; G D) vers yo quand e 
dans Sq. Cette solution y{x, e) est holomorphe en ses deux variables dans T) x Sq. 
En regardant un peu plus loin dans le developpement de y on arrive au theoreme 
suivant : 

Theoreme 11 La solution y du theoreme 10 adrnet la solution formelle y = 
Ylyn{x)e"' comme developpement asymptotique, uniforme pour x ^ T>, quand e 
tend vers 0. 

Nous avons deja montre que la vraie solution y{x, e) existe et qu'elle est bornee 
dans D. Par hypothese, on salt pouvoir determiner, pour tout N, les premiers 
coefficients de la serie formelle y : yoix), . . . ,yjv-i(x). Posons Risf{x,e) tel que 

y{x, e) = yoix) + eyi{x) ^ h £^~'^yN-i{x) + £^Rn{x, e) . 

On salt de la fonction Ri^{x,e) qu'il s'agit, pour tout e & Sq fixe, d'une fonction 
holomorphe et bornee en a; G D, et que 6^Rn{x, e) est bornee, uniformement pour 
X & D, quand e tend vers dans Sq. On souhaite montrer en plus que RN{x,e) 
est bornee sur D x 5*0, en utilisant une equation difFerentielle pour R. 

En remplagant dans (l) , on obtient une equation diflerentielle pour -Rjv, puisque 
tons les autres termes sont connus. Dans le cas = 1, pour lequel on sait en plus 
que Ri{x,e) est bornee sur D x Sq, 

e {yo + sRi)' = {x^f + eg) {yo + eRi) + h + e{yo + eRif P{x, e, ey + e^Ri) 

eR[ = {x^f + eg) R, + yog -y'o+ ~ °^ + y'P{x, e, ey) , 

(on ne notera plus les dependances en x et/ou e des fonctions la ou il n'y a pas 
d'ambigui'te) 

On obtient bien une equation differentielle en R-i{x,e), les fonctions y et yo en 
particulier pouvant etre considerees comme donnees du probleme. Cette equation 
differentielle est lineaire en Ri, mais on va se contenter de verifier que les fonctions 
coefficients verifient les proprietes du theoreme 10. 

Les fonctions f et g ne changent pas, leurs proprietes restent. La nouvelle 
fonction pour P, Pi, est nulle. Le nouvel h est 

, . , , h{x,e) — h{x,0) , 
hi{x, e) = yog -yo^ h y P{x, e, ey) , 
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une fonction holomorphe en x sur Dq, sauf eventuellement en a; = 0, a cause de la 
presence de h dans I'expression ; elle est aussi holomorphe en e dans 5*0, et reste 
bornee quand £ — > 0, d'apres ce qu'on salt de y et de la fraction en h. Enfin, hi 
est bien dans I'espace : 



hi{x,e) 



xPf 



< ho 



9 



xPf 



+ 



xPf 



\\yo\\ + 



sup 



dh{x, e) 



de 



\y\ 



P{x,e,ey) 



xPf{x) 



puisque toutes les normes citees existent bien. 

On pent du coup appliquer I'integralite du theoreme 10, la fonction Ri con- 
verge quand e tend vers 0, vers une fonction yi{x), uniformement pour a; € D : 
Ri{x,e) = yi{x) + eR2{x,e). La fonction Ri est, pour tout £ assez petit, bornee 
dans D. 

Rien n'empeche de refaire le meme raisonnement pour R2 a partir de I'egalite 
Ri = yi + R2- On a done, par recurrence directe. 



i<N-l 



< e^\\RN{x,e) 



ou Rn est, pour tout N, une fonction bornee en e uniformement sur D : y est un 
developpement asymptotique de y. 

Pour etre complet dans la demonstration, il faut remarquer que dans les hy- 
potheses H«, on demande que h{x,e) admette un developpement asymptotique 
en e quand e tend vers dans 5*0, ce qui n'est pas le cas a priori pour hi{x,e), 
puisqu'on n'a pas d'abord cette propriete pour y{x,e). Mais ce n'est pas un prob- 
leme : la demonstration du theoreme 10 demande seulement que I'on puisse ecrire 
h{x,e) = h{x,0) + ehs{x,e), avec bornee dans T) x 5*0, et cela, on sait deja 
pouvoir le faire pour y{x,e) — yo{x) + eRi{x,e). Ce qui suffit pour appliquer la 
recurrence. 



3.2.2 Au voisinage de 

Dans cette section, nous allons montrer qu'il existe des solutions pres du point 
tournant 0, tres exactement pour des x dans des domaines 

{ge'' < \x\ < p, arg(a;) G [^i, ^2]} , r = . 

Ces solutions seront aussi appelees solutions interieures. Comme il s'agit de pro- 
longer les solutions trouvees au paragraphe precedent, on conserve les hypotheses 
H«, et on fixe a nouveau (au moins dans un premier temps) un argument pour e, 
dans 5*0. 

Pour arriver au resultat de prolongement, nous regardons maintenant ce qui 
se passe au voisinage de a; = 0, en utilisant une hloupez. Ce changement de 
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variable donne une equation non singulierement perturbee permettant effective- 
ment d'analyser ce qui se passe au voisinage de 0. L'equation obtenue est appelee 
equation interieure. 



On pose 



Y ^ey ,et X ^ 



X 



sachant que e' = e^^^^^\ en prenant (par exemple) la valeur principale de la 
racine. En partant de £ e Sq, on arrive ainsi k s' E S'q. 
L'equation (l) transformee s'ecrit alors 

^ = (^X^fie'X) + e'gie'X, e)) Y + e'h{e'X, e) + j^Y'P{e'X, e, Y) . (ii) 

On va se placer dans un domaine Dq, qui contient une partie de la montagne 
qui prolonge celle de D vers et de ses deux vallees adjacentes, soit, a pen de 
choses pres, I'intersection d'un secteur ouvert centre en 0, d'amplitude et de 
taille p/e' et de la boule \X\ > g 
arbitrairement grande) : 



Xi\ (on se contentera d'une constante g 



-Dc^ <X eC / g<\X\< p/e' et aigXe 9,0 + 



A X donne, on definit le chemin Vx comme etant le chemin descendant le relief 
entre le point Xp/e' et X. Xp est pris dans I'ensemble des points «au-dessus » de 
e'X et de module p (cf. aussi fig. 8). 

On a defini le domaine Dc, de telle maniere que, pour tout X G Dc il existe 
un tel chemin Fx- 

Une hypothese necessaire pour continuer est que les trois limites suivantes ex- 
istent, au moins en tant que limites ponctuelles, pour tout X dans un domaine in- 

dependant de e incluant Vc, comme I'ensemble > \Xi\ et argX G 

et pour tout Y, \Y\ < S, pour un certain 5 fixe : 



Pe{X,e',Y) 



e'gie'X,e' 
e'h{e'X, e 



P{e'X,e'^+\Y) 



e'->0 



ip+l\ 



Go{X) 
Ho{X) 



Po{X,Y) 



Si on regarde ce que I'existence de telles limites signifie, on voit que ces hy- 
potheses sont equivalentes a : 
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Fig. 8: Domaine T)c 




1. g et h ont, au pire, un pole en x = 0. Si g a effectivement un pole, alors 
g{x, e) = go{e) + gi{e)x + g2{e)x + ... + gn{e)x^ + . . . + 



X 

en ecrivant g sous sa forme de serie de Laurent en ,t = (en mettant e 
en facteur la ou il est necessaire), dont tons les coefficients gk, k E Z sont 
encore des fonctions holomorphes pour e dans So, bornees quand £ — > ; 
d'ou, ensuite 

^ _ ^™£^ofi'-i(e) , lim^^o9-p-2{£) . lim^^o 5'-2p-3(e) , 
^o(A j + + -^^^ + . . . 

Si au contraire g est holomorphe en a; = 0, alors Go{X) = 0. On obtient les 
memes resultats pour h et Hq. II est clair que, dans tons les cas, e'g{e'X,e) 
converge uniformement vers Go{X)^ puisque tous les termes des series de 
Laurent restent bornes dans Dc : s'X et 1/X sont en effet bornes dans ce 
domaine. 
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2. P{x, e, ey) peut s'ecrire sous la forme x^Pi{x, £y)+eP2{x, £, ey), et on trouve 
dans ce cas-la que 

P,{X, e\ Y) = XPP,{e'X, Y) + e'P2{e'X, e, Y) , 

puis Po{X,Y) = XPP,{Q, Y). 



Pour commencer, nous allons montrer qu'il existe, pour X assez grand, une 
solution lo, bornee dans un tel domaine 'Dq, de I'equation obtenue en posant 
e' = dans (ll) 



dYo 
dX 



/(O) + G'o(X))yo + HoiX) + Y.^PoiX, Yo) (3.5) 



II s'agit la d'une equation irreguliere singuliere de rang p+1 (en I'infini), dont 
on sait qu'elle a des solutions asymptotiques. On redemontre ici sommairement 
ce resultat. 

D'apres leurs definitions, Go{X), Hq{X) et Pq{X , 5) / Xp+'^ tendent vers 
quand X tend vers I'infini. On peut done montrer qu'il existe une solution formelle 
en 1/X a cette equation(3.5) : ^o(-'^) = ~ /(o)xp — ' • Nous allons montrer, avec 
un point fixe, qu'il existe une et une seule solution a cette equation, qui tende vers 
en I'infini, dans certains secteurs d'amplitude inferieure a ; la serie formelle 
ci-dessus sera serie asymptotique a cette solution. On trouve d'abord que Yq est 
solution de I'equation 

Yo{X)^J exp(^^ («V(0) + G'o(«))ci«^ (^i/o(i) + >o(^)'Po(i,l^o(i)))c?i 

ce qui definit un operateur sur un sous-ensemble de I'espace des fonctions 
bornees dans D^o, qui tendent vers en I'infini. Les images de telles fonctions par 
cet operateur sont aussi dans cet espace. Puis, si Yq' et Fq" ^^^^ deux fonctions 
de cet espace on a 



, , , f f,^ i uPf(0)+Go(u)]du , , 

J^o(i^o^) - yo(Y,-) = / e^' \ V (r+ - y-) 

J Fx 



\\yo{Yo^) - yo{Y,-)\\ < 11^0+ - Y,-\\ ^ exp (^^(^j\u^ f{0) + Go{u))du) 



BP 

X (25 sup \Po{t,6')\+6snp -^{t,6') 
\s'\<s \s'\<s or 



dt 



61 



ce qui montre bien, puisque ^-{t,5)/t'P et Po{t,S)/tP sont bornees, qu'on a un 
operateur contractant pour tous les X accessibles, au moins pour I'ensemble des 
fonctions Yq majorees par un certain 6 : 

5 > 11^0^11= sup|yo^(x)|, 

\X\>s 

ce qui pour un () fixe est vrai pour tous les g assez grands, pusique les fonctions 
Yq tendent vers 0. De meme, si Hloll < ^5 alors ||3^o(^o)|| est inferieur a 6 pour un 
g assez grand, puisque Hq tend aussi vers en I'infini. 
Done lo existe bien dans un domaine de type "Doo- 

On remarque que si I'equation (ll) est effectivement non singulierement per- 
turbee, on souhaite montrer un resultat d'existence non trivial, puisqu'il s'agit 
de solutions dans un domaine non borne (parce que dependant de e), et avec des 
conditions initiales dependant de e (puisqu'on veut prolonger des solutions du 
paragraphe precedent). 

Nous allons demontrer I'existence de solutions Y de (ll) pour tous les X acces- 
sibles a partir d'un point Xp/e' fixe. Nous les cherchons sous la forme Y — Yg+e'Z. 
Alors 

dZ 



— ^[x^f{e'X)+e'g{e'X,e))z 



e'h{e'X,s)-Ho{X) f e'gie'X,e) - GojX) f{s'X)-m \ 
e' V ^' / 

, ^^, e'-'P{e'X,e,Yo + e'Z)-Po{X,Y,) , f^^^^ , ,^,^ P{e'X,e,Yo + e'Z) 
+ Yq + [^Yo^ + eZ j — 

(3.6) 

(on omet de rappeler systematiquement la dependance de Iq en X et de Z en X 
et e') 

On pose 



^^^^ ^ . e'g{e'X,e)-Go{X) ^ ^ J{e'X)-f{0) \ ^ e'h{e'X, e) - H,{X) 



e' £' / e' 

H{X, e', Z) ^yo^ ^^^^''^'^° + y^~^°^^'^°^ + {2Y,Z + e'Z')P,{X, e\ Y, + e' Z) 

+ P2{e'X, e, Yo{X) + e'Z) + {2Yo{X)Z + e'Z^) P,{X, e', Yo{X) + e'Z) 
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Avec ces notations, on prend 

Zh{X, e') = i {ey{x,, e) - Y,{x,/e')) + / e^^ ^'^^^ e') + H{t, e\ Z))dt , 
£ Jfx 

(3.7) 

ou a e' fixe, le premier terme correspond a une condition initiate aux limites du 
domaine Dc- 

On va utiliser un theoreme de point fixe avec des operateurs ressemblants a 
ceux du paragraphe precedent : 

7t '. Z) (Z Z, — > 0~C 

Z^H{Z), 

la fonction n{Z) etant telle que n{Z){X, e') = H{X, e', Z{X, e')) pour tout X e 
Dc, et 

Z : 0^ (Z 0^ — > Z, 
H^Zh 

la fonction ZH{X,e') etant definie par la relation (3.7). 

On pent noter que 'H{Z) contient ici encore quelques termes lineaires en Z, 
qui resteront toutefois bornes par rapport a ceux presents dans Fg. 

Les deux operateurs agissent a nouveau sur des sous-ensembles bornes d'es- 
paces IK et Z. Ces espaces fonctionnels sont les suivants : 



3-C = H[X) I H est holomorphe pour tout X e Ttc, et sup 

X 



H{X) 



XPf{0) 



est borne 



et 



Z — i^Z{X) / Z est holomorphe sur Dq, et tel que sup |^(-'^)| est borne | 
(les sup sont bien sur pris pour X dans le domaine Dc)- 

X 

On voudrait qu'avec les normes correspondantes W-Wj-^ et ||.||, I'inegalite suiv- 
ante soit vraie : 

\\H{Z,)-n{Z2)\\r^<CH\\Z,-Z2\\. 

Or, si \\Zk\\ < 5, 

\n{Z^){X,e') - n{Z,){X,e')\{X) < \\Z^-Z4 X {\Y^{X)\ + \e'\5) 



sup 

AS'\<5 



dY 



{X,e',Yo + e'S') 



{\Yo{X)\ + \e'd'\) + 2 sup \P,iX, e' , Yq + e'd' 

\S'\<S 
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Une condition suffisante est que, pour tout X dans Dc et tout e' assez petit la 
condition suivante est verifiee : 



sup 

\S'\<S 



dY 



{X,e\Yo + e'5') 



. (\Yo\ + \e'5'\) + 2 sup P,{X, e\ Fq + e'5') 

^ ' \8'\<8 

<M,|/(0)X^'| (3.8) 



Elle Test en particulier si 



sup 

|<5'|<5 



< M' 



!{e'X)X^ 



II suffit de revenir a la variable x pour s'apercevoir que cela correspond a une des 
hypotheses de croissance comparee du theoreme 10 : 



P{x,e,ey) Pi{x,sy) , P2{x,e,ey) 



f{x)xP 



+ e- 



xPf{x) 



est bornee . 



Cette hypothese reste vraie jusqu'en ,t = pour le premier terme puisque Pi{x, ey)/ f{ 
est holomorphe en ce point. En ce qui concerne le second terme, il reste borne 
pour tout X = e'C, C assez grand : e/x^ est bornee pour les x de cette forme, et 
P2/ f est elle aussi holomorphe jusqu'en x — Q. 
On a alors 



- n{Z2)\\<yi < {\\Yo\\ + \e'\6^M,\\Z, - Z^ 



(3.9) 



Par ailleurs. 



\Z{H,)-Z{H2)\{X,e')< [ |e^^(^)-^^(*)t^/(0)|citx ||ifi-i/2||j£; 

J Fx 

comme e'g{e'X,e) est une fonction bornee sur et que \ f{e't)tP\ = O (|t|^), ce 
qui implique Kl,\t\P < \Fl{t)\ < i^f-l^l^, on voit que 



\\Z{H,) - Z{H2)\\ < Cz.\\H, - H^Wj^ . 



(3.10) 



Or sup|jf|>^|Fo(-^)| + pent etre rendu aussi petit qu'on veut quand g croit et 
que e' est suffisamment petit. Comme par ailleurs Cz est une constante qu'on pent 
rendre independante de £ et X, et que est independante de e, les deux inegal- 
ites (3.9) et (3.10) permettent d'envisager que I'operateur ZoTi soit contractant, 
au moins pour les X dans un certain domaine Dc C {X / \X\ > g}. 

II reste a verifier que les images de ces operateurs (restreints sur des boules 
ouvertes de "K ou Z) sont bien incluses dans les espaces voulus. 

Si Z est une fonction de Z, dont la norme est majoree par un certain 6 > 0, 
la fonction image H{t,6', Z[t,6')) est bien dans I'espace "K : 
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► Ps est dans IK (puisque la condition (3.8) est verifiee). 

► Pi{e'X,Yo + s'Z) et l/XPP2{e'X,e,Yo + e'Z) sont holomorphes en toutes 
leurs variables, et done bornees pour X G Dc et \\Z\\ < S. 

► comme Pi est holomorphe en e'X, y compris au voisinage de 0, on pent 
ecrire P,{e'X, Yo + e'Z) = En>o(^'^)"^i,n(>o + e'Z), done 

^(Pi(£%yo + £'^)- A(o,yo)) 

= ^ [Pi,o{Yo + s'Z) - P,,o{Yo) + £'X^(£'X)-iPi,„(ro + e'Z) j 

\ n>l / 

= Z^(0, Fo) (l + oAl)) + X J2i^'X)--'P,,^{Yo + e'Z) 



n>l 



et puisque la somme dans cette expression est bornee par une constante 
C||Pi||, et que la derivee de Pi Test aussi, on pent majorer cette difference, 
dans laquelle intervient toutefois a nouveau un facteur X. 
Cela donne 

dPi 



\H{X,e',Z)\<\\Y,^-6 



dY 



(0,Vo) 



\XY,^\ 



C- \\Pi\ 



l-'o I 



IP 



2lW 



+ {2S ■ \\Yo\\ + e'S^) ■ \\P,{X, e', Yq + e'Z) ||^ . 

On voit que non seulement Ti.{Z) est dans I'espace !K, mais que sa norme pent 
etre choisie aussi petite que Ton veut, si on prend e' assez petit et g assez grand, 
et ceci pour tout 6 fixe. 

Dans I'autre sens, regardons quelle est Timage de la fonction nulle par Z. Cette 

fonction image est bornee si G{X,e') est dans !H et si —(^ey{xp,e) — Yo{xp/e')) 

est bornee quand e tend vers 0. 



G{X, e') = 



e'g{e'X,e)-Go{X) ^ ^^^ /(g^X) - /(0) \ ^ e'h{e'X,e) - HojX) 
e' e' J e' 



Or, d'apres ce qui a ete dit sur les consequences des hypotheses (H^q), 
e'g{e'X, e)^e'Y, gk{e){e'X)'^ + ^ 



fe>0 

e 



I ^,p-i g-2{e) ^ ^ ^, g-p-i{e) \ ^ g-p-2{e) ^ 



e 



X2 ) 

e'g{e'X,e)-Go{X) g_,{e) - g.,{0) 
^ 2^9k[e){eX) + — 



A:>0 



f ,p-i9-2 ^, 9-p-i \ 9-p-2{e)-c 



9-p-2{^) 



+ ... 
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Done , et de la meme maniere , sont des 

e' e' 
fonctions qui sont bornees dans Dq ; elles sont dans I'espace !K, et leur norme 

dans cet espace tend vers avec e. 

fle'X) — f(0) 

En revanche, la fraction n'est pas uniformement bornee pour 

e' 

X e Do- On peut cependant demontrer le lemme de maj oration suivant : 

Lemme 2 SoitF{x) une fonction holomorphe dans la boule ouverteT>F — B{0, p) . 
Alors, pour tout X tel que (s'X) & Dp, 

oil F est une fonction analytique ; de plus, si F est bornee, F est bornee elle aussi. 
On ecrit d'abord F{x) — ^„>o fn^"', puis on en tire 



n>l 



On a done F{x) = J2n>ifn^"' hien analytique en x, et qui est bornee 

pour X eVp si F I'est.'n 

fis'X) — f(0) 

Done la fraction peut etre d'ordre 0{X) ; mais parce que Yq{X) 



f{e'X)-f{0) 
e' 



tend vers quand X ^ oo au moins aussi vite que 1/X, la fonction Yq[X) 

f(e'X) — f(0) 

est, elle, bien bornee, et done X^YQiX) ; est une fonction de "K, elle 

e 

aussi. 

II ne reste done qu'une derniere condition a controler : on veut montrer que la 
condition initiale Z{xp/e' ,e') — ^^ey{xp,e) — Yo{xp/e')^ est aussi petite que Ton 
veut pour tout e' assez petit. Or y{xp, e) est borne en e, et YQ{xp/e') est equivalent 



Yo{xp/e') 



Ho{xp/e') ,pHo{xp/e') 



.0 /(O) {xp/e'r 



~ e 



sachant que Ho{xp/e') est au mieux un 0{e'). On est done assure que Z{xp/s',s') 
est d'ordre au pire e'^ quand e' tend vers 0, ce qui nous suffit largement. 

Done G{X, e') est une fonction de Di, on peut done ecrire 



\m\\ < \\yo\ 



s'g{s'X,s)-Go{X) 



XWo{X) 



fie'X) - /(O) 



+ 



e'h{e'X,e) - Ho{X) 
7' 



+ \Z{xpls! ,s! 
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Tous ces termes sont bornes, pour tout e' assez petit, et ils decroissent avec q — 

On choisit done 5 > 2||2^(0)||, puis on augmente la valeur de q et on fait 
deeroitre eelle de \e'\ jusqu'a ce que 

C^M^dlroll < 1 et que sup s', Z)||^ < bjCz. 

||Z||<25 

Dans ce cas, la composition des deux operateurs Ti. et Z est bien un operateur 
contractant. 

On pent alors enoncer le theoreme suivant. 

Theoreme 12 On suppose que pour I'equation (l) les hypotheses (H^q) et celles 
du theoreme 10, pour tout r > 0, sont verifiees dans un x-domaine accessible 
(avec le relief 'R^^), qu'on note d nouveau D C Dq. Soit Xi un reel assez grand, 
independant d'e si \e\ < Eq- Soit alors pour un tel e fixe des domaines 

Vc = \^X eC / g < \X\ < p/e' et argXe 

6i et 02 etant choisis, pour un e donne, tels que tout point X de T)c soit accessible 
(avec le relief ?R.(e'u^^^ f (0)) ) a partir d'au moins un point de la forme x/e', oil 
X eT) ; et 

= 2) u {x e C / (x/e') e Dc} . 

Alors il existe une solution holomorphe y de (l) sur D^i telle que ey est bornee 
quand e tend vers 0, et ey{Xe',e) = Y{X,e') tend, uniformement pour X dans 
T>c quand e tend vers 0, vers une fonction Yq[X) solution de I'equation (2.2). 

On salt qu'il existe une unique solution au probleme de Cauchy suivant : 
Y{X, e') est solution de (ll) et verifie Yixpje' , e') = ey{xp, e). On vient de montrer 
que cette solution existe sur la partie du domaine Dc accessible a partir du Xp/e' 
choisi et qu'elle reste dans un certain voisinage de Yq{X). Or, par unicite des 
solutions avec la condition initiale au point Xp, la solution de (l) y{x,e) pour 
X & T>,\x\ > p se prolonge avec la fonction ^Y(x/e',e) dans toute la partie de 
Dc accessible ; et comme tout Dc est accessible a partir des Xp/ e' (et meme 
d'un nombre fini de tels points), et qu'il ne pent y avoir qu'un seul prolongement 
dans ce domaine, y{x,e) se prolonge en fait avec -Y{x/e',e) pour I'ensemble des 
X e e"Dc. 

Ce qui clot la demonstration du theoreme. 

Cette solution y{x, e) est holomorphe en ses deux variables dans D^/ x Sq. Pour 
demontrer la dependance holomorphe en e', puis en e, de la fonction Y{X,e) cor- 
respondante, on precede comme pour le theoreme 10, en prenant un D'^ accessible 
pour un intervalle d'arguments de e' et en considerant des espaces de fonctions 
sur Dc X 5*0. On ne detaillera pas la demonstration du theoreme 12 modifiee ici. 



h,92 
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3.2.3 Developpement de la solution interieure 

La solution Y{X,e') trouvee tend vers Yo{X) quand £ — > 0, uniformement 
pour X G Dc- Admet-elle un developpement asymptotique de la forme Y — 

^„>Q F„(X)£:'", uniformement dans tout le domaine Dc ? 
La reponse n'est pas toujours positive. Si on tente de poser 

YiX,e')^ J2 Y^{X)e"' + e"'^'Y{X,e') 

0<n<N 

et de montrer que Y est bornee de la meme maniere qu'on a montre que Z 
etait bornee (theoreme 12), on se retrouve face a un probleme : an cours de la 
demonstration, on devait montrer, page 66, que G etait dans IH, ce qui n'est vrai 
que parce que Yq tend vers en I'infini. Mais on pent seulement montrer en general 
que Z est bornee dans Dc ; une demonstration par recurrence generale, valable 
sur tout le domaine TIq comme dans le cas de I'equation de Van der Pol, est done 
exclue. 

Par contre, on pent arriver a demontrer le corollaire suivant. 

Corollaire 13 La fonction Y{X,e') du theoreme 12 verifie la propriete suivante. 
Soit g]0, 1[, et g > un reel assez grand avec le domaine Dc qui en depend de 
la maniere habituelle . II existe une suite de fonctions Yq, Yi, . . . , Yn, . . . telle 
que, pour tout N , 

N-l 

YiX, e') - J2 yniX)e'" = e'^'z^iX, e') , 

n=0 

la fonction Zf^(X,e') etant definie pour X e "Dc, mais bornee par e'"^^"^^'^ quand 
e' tend vers et seulement pour les X tels que g < X < \£'\~'^ (on notera Dcfx 
I'ensemble des X bornes de cette maniere et qui sont aussi dans Dc). 

Nous demontrerons cette propriete par recurrence sur A^, en reutilisant une partie 
de la demonstration du theoreme precedent. On y traite de fait le cas = : 
Zi{X,£') est une fonction bornee, en particulier pour les X < l^'l^^. 

Mais auparavant. nous allons montrer certaines proprietes de majoration. 
On introduit dans ce but, pour tout m G Z la classe de fonctions suivante : 
F{X,s',y,yo) G si F s'ecrit comme une somme finie telle que 

F{X, e', y, yo) = J] Ei{X)Fi{e'X, e', y, yo) 
I 

ou 

► Fi{x,e' AjAjo) est une fonction bornee, holomorphe an voisinage de a; = 0, 
pour e' dans des secteurs Sq et pour \y\ < 5, |?/o| ^ 5; elle admet un 
developpement asymptotique uniforme quand — > 0. 
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► X "^Ei(X) est holomorphe et bornee independamment de e' dans Dcf^- 

Comme propriete de ces classes de fonctions, on pent noter que, si F e et 
G & J-'n, alors 

•k F + G est dans J^sup(m,n)- 

★ F X G est dans J-m+n- 

OF . .... 

•k —— est une fonction de J-jn, quitte a diminuer S. 

oY 

•k e'F G J-'m-i- 
Nous allons demontrer le lemme suivant : 

Lemme 3 Soit F{X,e',y,yo) une fonction de Tm- La fonction F definie comme 

F : {X,e ,y, yo) ^ 

e' 

est une fonction de la classe J-'m+i ■ 

La demonstration ressemble a celle du lemme 2. On ecrit Fi comme la somme 
convergente (pour s'X assez petit) 



Fiie'X, e', y, yo) = J](£'X)'=F^,,(£', y, yo) 



k>0 



On arrive alors a I'expression suivante (ou la reste convergente si e'X est assez 
proche de 0) : 

^fv ^' \ ^ ( VTP (r> \(^' v^*;-i_L -^^,o(^^ Ho) ~ ^1,0(0, y, yo) \ _ 
F[X,e ,y,yo) ^ 2_^EiX\ 2_^XFi^k{e ,y,yo){e X) H I; 

I \fe>i ^ / 

comme Fi o admet un developpement asymptotique uniforme en e' — 0, 

F{X, e'. y, yo) = E E^imHe'X, e' . y, yo) + E ■ 

\X-'^-^Ei{X)X\ et \X-'^-^Ei{X) \ etant bornees dans Dc^, la fonction F est done 
dans J-m+i- n 

Pour A'^ = 1, on reprend I'equation pour Z (3.6) vue plus haut, qui est verifiee 
par Zi : 

^ = f,{X)Z + G{X, e') + H{X, e\ Z) 
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avec 



/,(X) =f{e'X)XP + e'g{e'X,e) 
^^^^ ^ / e'gie'X,s)-Go{X) ^ ^^ J{e'X) - f{0) \ ^ e'h{e'X, e) - H,{X) 



y £' e' J e' 

H{X, e', Z) =Y^ PeiX.c'. lo + y) - P,iX. Yo) ^ ^^y^^ + e'Z')P,{X, e' , Y, + e' Z) 

P,{X,e',Yo + e'Z)-Po{X,Yo) ^ ^ y)ZP,{X, e' , Y, + e' Z) 

£ 

Pour simplifier au maximum le terme if, nous allons reorganiser I'equation, en 
utilisant le fait que Y est maintenant connue, comme fonction holomorphe bornee. 
On reprend la decomposition suivante de P^ donnee parmi les hypotheses : 

P,{X, e\ Y) = XPPi{e'X, Y) + e'P2{e'X, e, Y) , 

puis 

XPP,{e'X, Y) - Po{X, Yo) = X^ (Piie'X, Y) - P,{0, Y) + P,{0, Y) - P,{0, Y^)) 

= XP {P,{e'X,Y)-P^{0,Y)) 

+ XP{Y-Yo)J^ ^(Q,Yo + t{Y-Yo)yt- 

ce qui donne 

P,{X,e',Yo + e'Z)-Po{X,Yo) ^ ^^ P,{e'X,Y) - P,{Q,Y) ^ 

dPi 



+ X'zJ^ ^(o,Yo + t{Y-Yo)yt. 



On pent obtenir ainsi comme equation differentielle pour Zi 
dZi 



dX 

avec 



= h [X, Y{X, £'), Yo{X)) Z, + Gi (X, Y{X, £'), O) 



e\ y, yo) =f{e'X)X^ + e'g{e'X, e) + {Yo{X) + P,(X, e\ y) 
-rX^Y^iX) l'^(^0,yo + t{y-yo))dt 

G.iX,s',y,yo) = (^'^^^'^^f^^^^^ ^X^li^^^y^iX) 

e'M:'X.:)-H„{X) 
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C'est a partir de cette equation devenue lineaire en Z que nous ferons la demon- 
stration par recurrence. 

D'apres ce qui a ete vu plus liaut ces fonctions verifient les proprietes suivantes : 

► fi{X,e',y,yo) est de classe J-'p. 

► Gi{X,e' ,y,yo) est aussi de la meme classe J-'p. 

► Zi est une fonction bornee dans D^, done bornee dans Dcp,- 

On note Yi{X) la solution de I'equation de Zi avec e' = qui est bornee dans 
I^oo (^1 G -^o) ; on montre I'existence de cette solution comme on a montre celle 
de Yo, sachant que X'Pf^iX, 0, Yo{X), Yo{X)) et X-pGi{X, 0, Yo{X), 0) sont deux 
fonctions bornees dans Dcn- 

On pose ensuite 

Z^{X,e') = Y,{X)+e'Z2{X,e'). 

Nous aurons done besoin de completer le lemme 3, en regardant ce qu'on pent 
en deduire pour les classes Q^i^) comprenant les fonctions G{X,s',y,yo,Z), ou 
G est un polynome en Z dont les coefficients Gi de sont dans J-'k-im- 

L'interet principal de ces classes Q de fonctions est que, si on salt que X~"^Z{X, e') 

est une fonction bornee dans T>c^, alors G(^X,e',Y{X,e'),YQ{X), Z{X,e')^ croit 
an plus comme X^"*"""^ avec X dans ce domaine. 

Lemme 4 Soit G{X,e' ,y,yQ, Z) une fonction de Q^_2{Z). Alors on a les pro- 
prietes suivantes : 

1. G multiplie par une fonction de Ti est dans Q^^i^)- 

2. Le produit de G par Z est dans G^2 '^{Z). 
BG 

3. ^iX,e\y,yo,Z)eg'^-_T\Z). 

4- Une fonction G(X,e',y,yo,Zjn_i) de ^^_2(^m-i) pewi se transformer en 
posant Zm-i = Ym-i{X) +e'Zm- La fonction G{X,e',y,yo,Zm) obtenue est 
alors dans g^_i{Zm). 

1. Ce premier point est trivial. 

2. Pour demontrer cette propriete, on part de la definition : 

G{X, e', y, yo, Z) = J2 GiZ\ Gi{X, e', y, yo) G J'k-iim-2) 

I 

G{X, e', y, yo, Z)Z = ^ GiZ'+' = ^ Gi^iZ\ 

I I 

et Gi^i e J^k — l)(m— 2) J'k+m—2—l(m—2) 



71 



3. On precede comme ci-dessus 

dG 
dZ 



dG 

— {X, e', y, yo, Z)^Y. = + 

I I 



avec (/ + l)Gi+i e J^k-{i+i){m-2) = J^k-m+2-i{m-2) ■ D'ou la propriete annon- 
cee. 

4. La demonstration ici se fait en ecrivant Z^-i = i^m-i + £'Zm, ou Ym-i G 
On obtient alors un polynome en Z^ dont le terme en Z^ est une 
somme de termes de la forme 

oil Cj designe le nombre de combinaisons, et Gj est par definition dans 
J-'k-j{m-2) ; ce qui fait qu'on a en facteur de une fonction dans Tk-j{m-2)+(m-2){j-i)-i 
J-'k-i{m-i), pour tout /. On remarque en passant que G{X,e' ,y,yQ,Yjn_i{X)) 
pent se voir comme un element de J-'k- 

Ces proprietes servent essentiellement a demontrer le lemme suivant : 
Lemme 5 Soit G e Q^_2{Z) et H la fonction definie par 

H{X, e') = g(^X, e', Y{X, e'),YoiX), Z^-iiX, s')) ■ 

Alors il existe G £ Q^\{Zm) telle que 

I {h{X, e') - H{X, 0)) = e', Y{X, e'), Yo{X), Zm{X, e')) • 
On decompose 

G {X, s', Y, Yq, Zjn_i) — G (X, 0, F, Fq, -^m-l) 



s', y, Yo, Zjn) — 

+ 



£' 

G (X, 0, Y{X, e'), yo(X), Z^.Q - G (X, 0, 0), Y^{X), Z^_,) 

e' 



, G {X, 0, yp, yp, Zm-i{X, e')) - G {X, O, Fq, Fq, 0)) 



Dans cette somme, 

1. Le premier terme correspond, si on utilise le lemme 3, a un polynome en 
Zm-i dont le coefficient Hi[X,e' ,Y,Yq) pour Zm-i est dans Tk+i-i{m-2)- Si 
on remplace ensuite par {Y^_i + e' Z^)\ on retrouve le resultat de la 

propriete 4 ci-dessus. 

Le premier terme est done dans Q^\{Zm)- 
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2. On peut ecrire la deuxieme fraction de (3.11), a I'aide de la formule de Taylor 
avec reste integral 

YiX,e')-YiX,0) l^'^l^x, 0, [Yo + t{Y-Yo)] {X, s'),Yo, Z^.,{X, s'))dt . 

L'integrale st holomorphe en y et Iq; et c'est aussi un polynome en Z^_i 
qui reste dans Q^_2{Zm-\)- On le multiplie par 

e' 

= Yr{X) + e'Y2{X) + . . . + + e'^'-^Zm-i ; 

puisque e'-'^'^Yj^i{X) G JFq pour tout j, avec les points 1 et 2 du lemme, on 
salt que le produit est dans 

La deuxieme fraction est done (largement) dans 

3. La troisieme fraction se reecrit aussi avec une formule de Taylor 

Zm-l 

{X,e')-Zm-i{X,Q) 



e' 

^ dG{X, 0, y, Z) 
dZ 



(X, 0, Fo, lo, [Zm-l + t{Zm-l - Ym-l)\ {X, s')) dt . 



On integre un polynome en t dont les coefficients sont, d'apres la propriete 3 
du lemme precedent, des elements de ^^~™^^(-^m-i)7 on (avec la propriete 4) 
de G^-i^^iZm)- L'integrale est done dans G^-i^^{Zm)- On la multiplie par 
la fraction, egale a Z^. On se retrouve done dans Q^^^"^^^^ — Q^^.U 

On fait comme hypotheses de recurrence que pour m < n — 1, Z^, est solution 
de I'equation differentielle 

^ = /i (X, Y{X, £'), Y^{X)) Zm + Gm {X, s' , Y{X, £'), Yo{X), Zm-i{X, e')) 

oil Gm{X,e' ,y,yo, Z) est un element de Q^^~^{Zm-i)- La croissance de la fonc- 
tion Zm est majoree : X~"^+^Z^(X, e') est bornee dans Dc^- II existe alors pour 
tout m <n — 1 une unique fonction Y^ solution de 



dYr^ 



= h (X, 0, ro(^), Y^{X)yrn + Gm (x, 0, Yo{X), Y^{X), Yrn-l) 



dX 

telle que X-"'+'^Y^{X) est bornee {Y^ e T^-x). 
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Avec ces hypotheses, on pose 

Z^_,{X,e') = F„_i(X) +£X(X,£') . 

On remplace Z^-i par cette expression dans son equation differentielle (on ne 
notera plus les dependances de Y{X,s') et Yq{X)) : 



dX 



Gn-l{^, y, Yo, Zn-2) — Gn-l{^, 0, Yq, Yq, Yn-2) , s 
H - . {6. 11) 



On veut etudier la croissance de la fonction G„ on Ton pose 

n V V V \ - g^ Yq) - f\{X, 0, Iq, Yq) 

e' 

L'hypothese de recurrence est done que G„_i £ ^^l3^^(Z„_2). 

y , . MX,s',Y,Yo)-MX,0,Yo,Yo) 
bi on appiique le lemme 5 a ; , en considerant ji 

e' 

comme element de Qq{Zi), on voit que cette fraction est dans la classe (?^^2(^n-i) ; 
en multipliant encore par Z^-i, on se retrouve, d'apres le lemme 4 (propriete 2), 
dans g^t2~^{Zn-i). 

. V Gn~l[X,e' ,Y,Yo, Zn^2) — Gn-l[X,0,Yo,YQ,Yn-2) , , ^ 

Quant a , le lemme 5 nous 

e' 

indique qu'il s'agit d'une fonction de Q^'t^^l'^^ (Zn-i) , soit Q^t.2~^ {Zn-i) ■ 
Par consequent, G'„(X, e', F, Fq, ^n-i) e Qn-2~^i^n-i)- 

On applique maintenant la formule de variation de la constante pour I'equation 
differentielle (3.12) concernant Zn : 

Zn{X,s') = Z„ {x,e"'-\e') + / e^^W-^^WG„(t, s', F(t, e'), >^o(t), ^n-i(t, £'))t^^ 

(3.14) 

avec Fi{t) = /* fi (m, e', F(m, e'), Yo{u))du. Le chemin est lui deduit du chemin 
Fx, a I'interieur du domaine T>c^, descendant entre un point Xpe'^~^ (au lieu de 
Xpe'~ pour Fx) et X. 
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On en deduit que 



\Zn{X,e')\ < sup 



< sup |X|("-^)- 



/i(x,£',y(X,£0,>o(^,£')) 



(3.15) 



Cependant, cette demonstration doit encore etre completee : il reste a montrer 
que la ficondition initialez que I'on choisit egale a 

Z„ {x/^-\e') = ^ {ey{x/^, e) - Fol V') - ■ ■ ■ - V„_i(a;,£"^-^)) 

pour assurer la continuite entre solution interieure et exterieure n'est pas plus 
grande que les autres termes de I'inegalite (3.15). Cette propriete sera demontree 
au paragraphe suivant (coroUaire 16) ; pour le montrer il faudra chercher a etendre 
le domaine de la solution exterieure. 

On remarque que I'on construit avec cette demonstration du coroUaire une 
solution formelle Y[X^e') pour I'equation (ll), qui a une propriete ressemblant 
a celle d'une serie asymptotique pour la vraie solution Y{X,e'). Mais si on ne 
controle leur croissance de maniere interessante que dans Dc^i? tous les coefficients 
Yn[X) existent dans le meme domaine que ^0(^^)5 c'est-a-dire pour X e T)^. 

Comme I'equation (ll) n'est pas singulierement perturbee, on peut encore com- 
pleter le theoreme 12 et son premier coroUaire avec le coroUaire suivant : 

CoroUaire 14 (Principe de prolongement analytique) En plus de V ensem- 
ble des hypotheses du theoreme 12, supposons que la fonction Yq existe non seule- 
ment sur Dq, mais peut etre prolongee dans un domaine, accessible ou non, 
'Dp G Do, oil Dp est un domaine borne, tel que Dp fl D^ 7^ 0. On suppose 
aussi que toutes les fonctions g, h et P ont un prolongement analytique dans ce 
domaine. II existe alors une solution Z{X,e') holomorphe sur Dp qui est un pro- 
longement de la solution Y{X,e') du theoreme 12. Cette solution admetY comme 
developpement asymptotique uniforme sur Dp. Dans le domaine Dp U Dq, cette 
solution prolongee Y{X,e') tend uniformement vers Yq{X) quand e' tend vers 
et verifie encore le coroUaire 13. 

En particulier, si E Dp, on peut prolonger les solutions correspondantes de 
(l), y{x,e) = ^Y{x/6',6), jusqu'a un voisinage dex = 0. Ce voisinage, cependant, 
n'est que de taille autour de 0. 

Par ailleurs si on est sur que y n'est pas exponentiellement grande (em) dans 
ce voisinage, elle est a priori de taille 1/e. 
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La demonstration de ce coroUaire ne pose pas de probleme. 

Soit Xq un point de DpflDc- D'apres le theoreme, on a une solution holomor- 
phe Y{X, s) sur Dc, prolongeant y{x, s), et en particulier, Y{Xo, s') est bien defini. 
L'equation (ll), regulierement perturbee, a une solution holomorphe sur Dp, avec 
la condition initiale Z{XQ,e') — Y{XQ,e'). Cette solution admet un developpe- 
ment asymptotique ZnS'^, qui est necessairement egal a la serie formelle Y, par 
construction : les K„ et Z„ sont definies par recurrence comme solutions d'equa- 
tions differentielles (les memes pour Yn et quel que soit n) et ces fonctions 
admettent la meme condition initiale en Xq, done Y^ = Zj^. Par unicite, Z est 
bien entendu un prolongement de Y. 



3.2.4 Connexion des developpements asymptotiques 

On conserve dans ce paragraphe toutes les hypotheses des paragraphes prece- 
dents. 

En regardant alors le resultat du theoreme 11 et du corollaire 13, on s'apergoit 
qu'on n'a pas, pour I'instant, de moyen de donner une bonne estimation de la 
solution y{x,e) trouvee pour certains x au voisinage de : |£'|''^~'* < < \xp\ . 
Or on aimerait avoir ime telle estimation partout ou la solution est definie. 

Nous allons voir que la solution formelle exterieure reste un developpement 
asymptotique de la solution pour ces x, en adaptant les demonstrations du para- 
graphe 3.2.1. Les hypotheses H^q nous indiquent que si|x|>|£'|'',0<i/<l, 

\gix,e)\ < Bg\e'\~'' 

\h{x,e)\<Bh\s'r 

puis, comme / est bornee en et que h{x, 0) pent etre d'ordre 1/x en x — 0, 

h{x, 0) 



\yo[x) 



<By\ 



xPf{x) 

On reprend les memes definitions pour les espaces et les normes associees, et pour 
les operateurs entre ces espaces, qu'a la page 55, sauf que tout est cette fois defini 
pour X dans le domaine 

d;: = |x e Do / ^ = X eT>c et K\xp\ > \x\ > -^kT} 

Dans un tel domaine, les majorations pour les fonctions derivees seront les suiv- 
antes (quelles que soient les normes choisies sur des espaces de fonctions holomor- 
phes) : 

► si \f{x)\ < A pour tout x e D'^, alors dans tout domaine [x < k), et 
pour un certain B > A, < -j—^. II ne s'agit la que d'un cas particulier 

de la formule de Cauchy. 
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► si pour tout i^, e'^\h{x,e)\ < A, pour tout x e et e assez petit, alors on 



peut majorer la difference h{x,e) — h{x,0) : 



^h{x, e) — h{x, 0) 



< 



A 



En effet, loin de x = 0, h est holomorphe en x et admet un developpement 
asymptotique en e uniforme en x, ce qui fait que {h{x, e) — h{x, 0)) est borne 
par Ce ; au voisinage de x = 0, on h peut avoir un pole, on ne considere que 
la partie polaire de son developpement, qui nous empeche d'avoir le meme 
resultat : 



'polaire 



{x,e) 



h. 



X 



H h fh-nx"" + e 



h. 



n-l 



X 



H h 



h 



X 



+ 



Ici n est tel que ni' < m{p + 1) pour tout v < 1, done n = m{p + 1), et 
doit etre tel que Nu < {p + 1) +m{p+l), done N — n+p+1 ; en continuant, 
on peut avoir dans hpoiaire des termes en £^x~^^~^~", etc .... En calculant 
alors la difference h{x,e) — h{x,0), on trouve alors le resultat, qui est vrai 
pour toute fonction h verifiant la condition. 

Note : Dans tout le reste de ce paragraphe, on omettra de noter les constantes 
K > 1 successives decroissantes qui interviendront, et les constantes multiplicatives 
A, B, etc. 

On obtient comme majorations 

g{x,e) 



\\h\\^ < 



sup 



xPf{x) 



< e 



mW < \^\ 



et aussi 



< 



-pv 



\-1v 



alors que 



|P||^ reste bornee quand e tend vers 0. 



Cela va nous permettre de montrer le theoreme suivant, qui prolonge le theoreme 
11 : 

Theoreme 15 Si les hypotheses H^q et celles du theoreme 11 sont verifiees, alors 
pour tout V g]0, 1[, la solution y{x,e) et la solution formelle y sont reliees dans 
le domaine (defini ci-dessus) par la relation : 



N 



y[x,e} 



n=0 



X]E' 
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Nous commencerons par demontrer le cas = 0, en reprenant effectivement les 
memes operateurs que ceux qu'on trouve en page 55. 



Alors, si \\zk\\ < \s\ 



-2u 



{s Zk continue heureusement a rester aussi petit que 



necessaire), on obtient d'abord I'inegalite suivante, deduite de (3.1) : 

\\n{zi)-n{z2)\\ < |£| • ii^i-^all X 

{hh + 2(lbo|| + \e\\\zk\\) ■ \\P\\sj + \s\i\e\U\\ + ||yo||)'C||P||^)) 

< |£| [\er + i\er + kr"'i + m^r + kr"'i'] i^i - -^211 



\\'H{z,)-'H{z,)\\<{\s 



+ \e 



2-2u 



+ £ 



4-4v 



et on garde pour I'autre operateur I'inegalite (3.2), independante de e, vraie quelles 
que soient les normes des 

\\Z{H^)-Z{H2)\\<M.\\H^-H2\\^ 

ce qui permet d'envisager que la combinaison des deux operateurs soit contrac- 



tante, si les espaces d'arrivee sont bien les bons. Or, si \\z\\ < \e\ 
(3.3) 



^2u J) 



d'apres 



\z\ 



I ■ 1^1 

|l-3i/ 



-2u _ 

+ k 



gWsj + \\yo + £z'^ 
■ k + £ 



\\n^)\\sj < k 

< \e 

< \e 

< SUp(|£|-M£|^-=^'^) 

et dans I'autre sens, si ||-f^||j^ < sup^jsj"'^, |£|^~^'^^, on reprend I'inegalite (3.4) 

h{x, e) — h{x, 0) 



\Z{H)\\ < M 



yo9 - 2/0 + 



< 

< \£ 



( I ^\—^ I ^ I 
I |£| |£| 



+ M\\H\\^ 



+ k 



-2v\ 



+ k 



+ k 



-2u 



Ce qui permet d'appliquer effectivement le theoreme du point fixe et de demontrer 
le theoreme 15 pour le cas = 0. 

Pour le cas general, la methode est exactement la meme que pour le theoreme 
11; on obtient les memes equations pour les la seule difference est que les 
fonctions ne sont plus bornees independamment de £, mais sont bornees par une 



78 



puissance de e. Ainsi, en reprenant les notations de la demonstration, 



xPf 



< 



\yo\\ ■ hW^ + 



< \£\ \£\ 



-2u 



XPf 

-2u I 



+ 



sup 



h{x, e) — h{x, 0) 



+ NMl^lli3 



\e\ sup I \e 



1 



+ k 



-2u 



\\hi\\^ < \e\ 

Comme Ri verifie I'equation 

eR[ = {xPf + eg)Ri + hi 

on voit que ||i?i|| < et ||yi|| < puis avec £i?2 — Ri — Hi, 

\\R2{x,e)\\ < = \\h,\\^ + ll^lli^llyill < 

Par recurrence on arrive alors sans difficulte au resultat complet du theoreme, 
puisqu'on aura toujours 

hn-iix,e) -hn-i{x,0) 
hn{x,e) = g{x,s)yn-i{x) 



et 



|-Rn(2;,£)|| < ll^nll^^, comme |||/„(x)|| ^ ■ 



Ce theoreme de prolongement de la solution exterieure va nous permettre de 
finir la demonstration du paragraphe precedent, en assurant la connexion a tons 
les niveaux n entre la serie formelle interieiue Yl^n{x/s')e'^ et la serie formelle 

exterieure ^?/„(,x)£". Nous appellerons cela le 

Corollaire 16 Soit N un entier positif et n E ]0, 1[. Si les N premiers termes de 
la serie {Yo{X),Yi{X), . . . ,Yisf-i{X)) construite au corollaire 13 existent et sent 
majorees jusqu'd I'ordre N — 1 comme le prevoit ce corollaire, alors on a encore 



N 



n=0 



ouYN{X) = lim,,^oZN{X,e'). 

II est clair que ce theoreme suffit a completer la demonstration par recurrence du 
corollaire 13, puisque la condition initiale pour Zjsr+i (qui correspond au premier 
membre divise par s'^"*"^) a alors au plus la meme taille que celle prevue pour la 
fonction elle-meme. 

Soit z/ < /i, avec pourtant u peu different de /i. II est possible d'ecrire pour 
tout n > 0, d'apres le theoreme 15 (on rappelle que e'^~^^ — e). 

Y (xpe'^''^,e'^ = ey {xps'^, e) 



79 



k=0 



On decomposera y jusqu'a un uq tel que (no + 1) — (no + 2)/i > {N + 1) — N{1 — /i) . 
On connait une majoration des sur Tli, : 

Done les yk{x), qui sont independants de e, ont en a; = au plus un pole d'ordre 
{k + l)ip + 1). Autrement dit, Uk (pour tout k < n) qui est meromorphe en 
comme toutes les fonctions de I'equation (l), s'ecrit aussi, au voisinage de 

yk{x) = + yk,-^k+i)(p+i)+ix-^'''^^^^''^^^+'' + ...+ yk,o + ■■■ 

L'hypothese (de recurrence) principale du corollaire est 

7V-1 

qui implique 

k=0 1=0 

Le premier membre de I'inegalite ci-dessus correspond aussi a 

par definition de Z^, qui est majoree dans Dcad d'apres Fhypothese de recurrence. 
On pent en deduire une definition de Yn (xps"'"^) a un o (^|£'|('^~i)(^~i)) pres. En 
X , formellement, on a d'abord 

Y^{X)=\im^(Y{X,e')- Yl e"YkiX) \ . 

^ ^ ^ \ ik<Ar-i / 

Puis, 



no 



N-1 



/|iV|./|(M-l)(7V-l) 



< \£ \£ 



e'^YN{X) = Terme d'ordre e'^ dans ^s'^+^ykiXs') 

k>0 

= Terme d'ordre e'^ dans ^ 



•fe+i 



fe>0 



l>0 



ce qui mene a 

Yn{X) ~ 5^y,,iv-(,+i)(p+i)X^-('=+^)^+^) : 



fc>0 
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comme Yj^ est solution d'une equation differentielle avec une irregularite singuliere 
en I'infini, on salt que ce developpement formel est effectivement un developpe- 
ment asymptotique pour Yn{X). 
D'ou 

k<no 

On pent prndre no arbitrairement grand, done dans la difference 

no N 



k=0 1=0 



il ne reste, a un o (|£'|^"^^|£'|^'* ^^^) pres, qu'une somme finie de sommes conver- 
gentes : 



_fc+i 



yk,n {xpE^'y 



k<no \n>N-{k+l)(p+l) 

Or tons les termes de cette somme sont d'ordre, en e', au moins 

(k + l){p + 1) + (^N + 1 - (k + l){p + 1)^1^ 

^{k + l)(p +l){l-ii) + {N + l) + {N + - 1) 
^(N + 1) - N(l - + i^)(kp + p + k). 

Le dernier terme de cette expression etant strictement positif, I'ordre est superieur 
a (A^ + 1) - A^(l - Ce qui s'ecrit 



N 



k=0 



ce qu'on voulait demontrer. 
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3.3 Existence de solutions bornees jusqu'en x = 

On reprend I'equation originale, en y introduisant explicitement un multi- 
parametre : 

a{£) = (^ao{e), ... , ai{e), ... , ap-i{ej^ 
qu'on supposera avoir une serie asymptotique 

oo 
n=0 

On notera aussi 

p-i 

a{x,e) = 'y^^x^ai{s). 

1=0 

On considere done un eas partieulier de I'equation (l), qu'on ecrira par exemple 
sous la forme 

ey' = [x^ f{x) + eg{x., e, a))y + £, a) + ey'^P{x, e, a, ey) (ill) 

on encore 

ey' = {x^fix) + eg{x, £))y + h{x, e) + ey^P{x, £, ey) + Q{x, e, a, y) , (ill') 
ou on suppose que Q comporte les termes suivants : 

Q{x, e, a, y) — a{x, e)Ro{x) + eRi{x, e, a) + eyR2{x, e, a) 

+ ey'^R^lx, e, a, ey) + x^R^lx, e, a) , (3.16) 

avec i?o(0) ^ 0. Les deux premiers termes sont issus de h, qui doit done avoir une 
forme particuliere ; les termes i?2 et Rs viennent de 51 et P respectivement, dans 
le cas le plus general. 

Comme on va s'interesser maintenant a des solutions y{x, e) qui soient bornees 
sur un voisinage entier de 0, il parait necessaire que cette fois toutes les fonctions 
qui interviennent dans I'equation soient holomorphes en x — 0. 

Nous allons dans cette partie tenter de montrer que, dans certaines conditions, 
il existe de vraies solutions d'une equation differentielle qui existent jusque dans 
un voisinage d'un point tournant. Pour cela, nous allons montrer que les vraies 
solutions qui existent sur des montagnes (cf. 3.2.1) peuvent parfois etre prolongees 
jusqu'au point tournant, en se servant des resultats du paragraphe 3.2.2. Ensuite, 
en jouant sur un (multi-)parametre, on fera en sorte que ces differents prolonge- 
ments aient la mcme valeur en ,t = ; cela signifiera que tons ces prolongements 
correspondent a une seule solution de I'equation, qui existera done sur la reunion 
des domaines d'existence de ces solutions. 
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3.3.1 Existence d'une solution formelle continue en 



Dans le cas ou I'equation (lir) depend lineairement de a{x,e), elle pent se 
reecrire 

ey' = {x^f{x) + eg{x, e))y + h{x, e) + ey^P{x, e, ey) + a{x, e)Q{x, e, y) (3.17) 

Dans ce cas precis, on pent montrer qu'il existe une solution formelle a I'equation 
qui est continue en 0, a condition que (5(0, 0, 0) soit non nul. Nous allons montrer 
que I'equation (3.17) possede une unique solution formelle (^, a) qui s'ecrit comme 
serie formelle 

p-i 

y = ^""vM et a = ^ e^-anix), a„(a;) = ^ ak,nx'' , 

n>0 n>0 k=0 

OU les yn eux-memes sont des fonctions de x qui sont soit holomorphes en x, soit 
au moins admettent, quand x tend vers 0, une serie asymptotique bornee en 0, 
notee Tons les a„(a;) sont, quant a eux, des polynomes de degre inferieur ou 
egal a p — 1 en x. Pour obtenir ce resultat, il suffit essentiellement de remplacer 
dans (3.17) y et les ak par leur serie formelle, les fonctions holomorphes en e : g, 
h et P par leur serie, puis d'identifier les coefficients de 
On ecrira la fonction Q sous la forme 



Q{x,e,y) = ^Qn{x,y)s'' , 



n=0 

et on decomposera de la meme maniere la fonction P{x, e, ey) par rapport a e. 
Ce type de notation vaudra aussi pour h{x,e) et g{x,s). 

On commence pour n — par : 

= xPf{x)yo{x) + ho{x) + ao{x)Qo{x, y) 

(on remarque que Qo^x^y) = Q{x,0,0)). D'ou 

ho{x) + ao{x)Q{x,0,0) 
yoi^) - (3.18) 

On souhaite que yo soit continue en 0. II faut et il suffit pour cela que la valuation 
de hQ{x) + aQ{x)Q{x, 0, 0) soit superieure ou egale a p — 1 . Ce sera effectivement 
le cas si (et seulement si) on choisit pour ao(^) Is polynome obtenu par division 
des series en x suivantes : —ho{x) par Q{x, 0, 0). Cela nous donne bien une et une 
seule valeur possible pour qq puis yo- 
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Supposons a present que Ton ait montre I'existence (et Tunicite) des n premiers 
termes des series y et a, done regarde les coefficients de pour k <n — l. Alors, 
on obtient I'egalite suivante pour le coefficient de : 

nyn-i{x) = x^f{x)yn{x) + ^ gk{x)yn-i-k{x) + hn{x) + Pn-i{x,y) 

k<n-l 

+ ^ak{x)Qn-k{x,y) 

k<n 

On pent verifier que Pn-i{x, y), le coefficient de e'"~^ dans y'^P{x, e, ey) ne depend 
que de x et de yQ,yi,... ,yn-i- De meme, Qk{x,y), d'apres la forme de Q, ne 
depend que de a; et de |/o, . . . , yk-i- 
Cela permet de trouver 

nyn-i - ^ gk{x)yn-i-k + Kix) + Pn-i{x, u) + ^ ak{x)Qn-k{x, y) 

k<n—\ k<n 

xPf{x) 

Or I'expression 

nyn-i- 2^ gk{x)yn-i-k + hnix) + Pn-i{x,y) + ^ ak{x)Qn-k{x-,y) 

k<n—l k<n—l 

pent s'ecrire (au moins formellement) comme une serie en x si on remplace tons 
les yk pour k < n par leur serie (formelle) respective. On choisit alors le polynome 
an{x) tel que les p premiers termes du numerateur de la fraction ci-dessus soient 
nuls, i.e. que y„ reste bornee en x = 0, a nouveau en effectuant une division de 
series. 

II existe un unique polynome a„ de degre inferieur ou egal a p — 1 qui verifie 
la propriete voulue, puis le ?/„ calcule est bien sur unique lui aussi. 

Nous pouvons ainsi construire des series uniques a et y, solutions formelles de 
(3.17), telles que tons les coefficients de y restent bornes quand x tend vers 0. En 
revanche, rien ne prouve I'existence pour tout e assez petit d'une vraie solution 
y{x,e),a{x,e) a I'equation differentielle qui soit continue et bornee quand £ — > 
dans un voisinage de 0. 

Pour le cas non lineaire en a de I'equation (ill'), la methode de construction 
est analogue, mais plus complexe. Cependant, dans le cas qui nous interesse ou 
Q peut s'ecrire sous la forme (3.16), I'existence et I'unicite de la solution formelle 
{y,a) reste assuree (pour une demonstration complete, voir [4]). 
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3.3.2 Existence de vraies solutions holomorphes jusqu'en 

Theoreme 17 On suppose que pour I'equation (ill), en plus des hypotheses du 
theoreme 12, on a que les deux fonctions g et h sont holomorphes en 0. Alors 
il existe de vraies solutions (y,a), oil a s'ecrit necessairement a{x,e) = ao{x) + 
ai(x)£ + 0(e'^) (les deux termes Oq et a\ etant imposes), telles que les y{x,e) 
definies, bornees sur une montagne verifient en plus les proprietes suivantes : 

1. 1/(0, e) est definie et bornee quand £ tend vers 0. 

2. De plus 1/(0, e) - i/(0, 0) = 0{e) (ou o{e) ). 

D'apres le theoreme 12, il existe des solutions y{x,s) de I'equation (ill) dans 
des domaines tels que x > \Xi\e' , ce qui ne permet pas encore de les definir en 0. 
Pour arriver jusqu'en 0, vu le coroUaire 14, il faut encore que Yo(X) puisse etre 
prolongee holomorphiquement jusqu'en X = 0. 

Or precisement, dans le cas qui nous interesse, ou g et h sont holomorphes en 
0, on a vu que Gq{X) = et Ho{X) = 0. L'equation pour Yq est alors partic- 
ulierement simple : 

(3.5) ^^{x^m))Yo + Y,'Po{X,Yo), 

et la solution Yq qui s'impose est une solution evidente : la fonction nulle, la seule 
solution de cette equation qui tend vers en I'infini dans un tres large secteur ; 
celle-ci existe bien sur en X = 0, et on pent utiliser le corollaire 14 qui assure 
alors I'existence de y{x,e) jusqu'en a; = 0. 

Mais ce resultat ne nous suffit pas, la solution y{x,e) obtenue n'etant pas 
bornee quand £ — > 0. Pour aller plus loin, il faut regarder exactement les premiers 
termes de la solution correspondante Y{X,e'), que dans les limites du corollaire 
13 on mettra sous la forme 

oo 

Y{X,e')r.J2^n{X)s'\ 

n=0 



Montrons que les p + 1 premiers K„ sont tons identiquement nuls. 
On salt qu'on a une solution yo{x) pour e = qui est continue en x = 0. 
L'equation pour y^ nous donne 

h(x,0,do) 

yo{x) = — , . 

j{x)xP 

done h{x^ 0, oq) est divisible par x'^. Si on pose a — aQ-\- ea, la fonction h pent 
se mettre sous la forme h{x,e,a) — xPhi{x) + eh2{x,e,a). Dans toute la suite, 
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on supposera qu'on a deja ecrit a sous cette forme et transforme g, h et P en 
consequence pour obtenir effectivement I'equation (ill) : 

ey' = (x^f{x) + £g{x, e, ea))y + x^hi{x) + eh2{x, s, a) 

+ ey^ (x^Pi{x, ey) + eP2{x, s, a, sy)j (ill) 

Le resultat pour I'equation en Y, (ii), est le suivant : 

^^(^X^f{e'X)+e'g{e'X,e,d.))Y 

+ sXPhi{e'X) + e'sh2(e'X, e, a) + Y'^Pe(X, s', a, Y) . (3.19) 

Done, par recurrence finie, en faisant I'hypothese de recurrence Yn-i{X) = 0, 
pour Yn, n allant de a p, les equations sont equivalentes a 



dYn 
dX 



= (xv(o))rn. 



La fonction nulle est la solution de chacune de ces equations qui verifie la condition 
initiale voulue ; c'est-a-dire, quand e' tend vers 0, 



e 



quand e tend vers 0, —y{xpe'^,e) tend vers et les autres termes a droite sont 
nuls par hypothese de recurrence, d'ou la condition 

lim Yn{X) = 0, Vn <p. 

X— >oo 

Pour n = p + 1, on obtient en revanche 

rIY 

v+i _ ixPf{0))Y,+, + XPh{0) 



dX 

dont la solution qui nous interesse verifie 



X 

c'est-a-dire 



lim Yp+,{X) = lim yoixs'") + o (s"-'') = 2/0 (0) 

f— »0O £'— »0 \ / 



Vi(^) 



■/(o) 
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On en conclut que Y{X, e') est egal a 

e') = eY{X, e') ~ + e%+^{X) + ...), 

d'ou Ton deduit que 

y{e'X,e) = ^^^^ = Y{X,e') 

£ 

reste bornee en x = (et dans un petit voisinage de taille e' autour de ce point) 
quand £ — > ; ce qui demontre le premier point du theoreme. 

II reste a montrer que I'equivalent de ?/(0, e) — ?/(0, 0) est d'ordre e au moins. 
On continue dans ce but le calcul des premiers termes de la serie formelle de Y . 

Revenons d'abord a I'equation (ill). Comme il existe une solution formelle, on 
sait calculer yi{x) avec I'equation suivante : 

y'^{x) = xPf{x)y,{x) + hix, 0, ai{x)) + yl{x)x^P^{x, 0) + yo{x)g{x, 0, 0) , 
done pour que yi soit continue en 0, il faut et il suffit que 

y'oi^) ~ h2{x, 0, ai) — yo{x)g{x, 0, 0) soit de valuation p au moins. (3.20) 
Pour Y, on trouve comme equation differentielle non singulierement perturbee 

^ = (^XPf{e'X) + e'g{e'X, e, ea)^Y + XPhi{e'X) + e'h2{e'X, £, a) 

+ eY^P,{X,e\S,,Y) 

On ne cherche que les termes Yq, Yi, ... Yp, done dans I'equation ci-dessus, on 
tronquera toutes les series pour ne garder que les termes de degre en e' strictement 
inferieur a p + 1. En particulier, le terme eY^P^{X,e',a,Y) pent etre oublie, et 
on obtient pour ces fonctions Yj des equations difFerentielles lineaires. 

Si on cherche a identifier les deux series formelles, celle pour y et celle pour Y, 
on s'apergoit que les Yj (pour j < p) ne devraient dependre que de yo, et etre de la 
forme Yj = \jXK Nous allons demontrer que les equations pour les Yj admettent 
effectivement de telles solutions. A partir de I'equation pour Y , 

^^jX^-'e'^ = [x^fie'X) + e'g{e'X, 0, 0)) ^ \,X^e'^ 

+ XPhi(s'X) + s'h2{e'X, 0, ai) [mod s'^+^J , 
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en ne regardant dans I'egalite ci-dessus que les termes dont le degre en e' est 
inferieur ou egal a p. 

V p 

^ e' (J2 ^jjX^'^s'^ - h2{e'X, 0, oi) - g{e'X, 0, 0) J] XjXh'^^ 

j=l 3=1 

p 

= (^f{e'X) J2 + H^'X)) [mod e'"^^] . 

j=i 

A gauche, on ne depasse pas le degre p — 1 en X, alors qu'a droite, on commence 
an degre p ; la seule possibilite d'avoir cette egalite est done que les deux membres 
soient nuls. 

A droite, cela donne (toujours en series tronquees) 

E ^ -^^^ = yo(s'X) [mod e''^'\ 

puis a gauche, on retrouve 

e' (y'^ie'X) - h2{e'X, 0, a{) - g{e'X, 0, 0)yo{e'X)^ 
qui d'apres la propriete (3.20) est de valuation p en e'X, done de valuation p + 1 



en e'. 



Ces Yj doivent verifier en plus une condition initiale, qui est, quand e' tend vers 
0, 



1 ^"^ / 



£'• 

k=0 



soit, si ^0(2;) = 'E.^^kx'', 
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Les solutions Yj trouvees formellement sont done bien aussi effectivement les so- 
lutions des equations differentielles eorrespondantes avee les eonditions initiales 
voulues. 

Pour y{x,s), on en deduit alors 

y{e'X, e)=Xo + eXiX + ■■■ + e"'\pX^ + 0{e) 
y{Q, e) - y{Q, 0) ~ £^^+1(0) voire plus petit, si Fp+i(0) = 0. (3.21) 

3.3.3 Une estimation preliminaire 

On etudiera dans la suite de cette partie I'equation (ill') dans laquelle on a 
deja ecrit a = ao + ea. 

ey' = {xPf{x) + eg{x, e))y + h{x, e) + ey'^P{x, e, ey) + Q{x, e, ea, y) 

la fonction Q etant holomorphe en toutes ses variables ; Q s'ecrit comme dans 
I'egalite (3.16). 

On note encore Q\ = i?o(0) 7^ 0. 

On notera 7^ un chemin descendant le relief a partir d'un sommet de la Z— ieme 
montagne jusqu'en x — (on rappelle qu'il nait p+1 montagnes autour du point 
0). 

Nous aurons besoin plus loin de I'estimation suivante : 
Lemme 6 Quand e tend vers 0, 



2m 
p+1 



/exp ( - I ff{t)dt + r{x) \ ^^(x, e, ea, y)dx ~ £Ck£'^~^^ exp [l{k + 1) 
oil Ck ne depend pas de I. 

II ne s'agit dans ce lemme que d'lm variante de la methode du point col. La valeur 
de I'integrale est essentiellement concentree la ou jjP f{t)dt est nuUe, c'est-a-dire 
en a; = 0. Comme exp(r(a;)) et f{x) sont holomorphes et non nulles au voisinage 
de X = 0, I'integrale est peu differente de 



( 1 /(O) \dQ ^ 

exp — x^^ - — (x,£,sa,y)ax 

e p+lj d&k^ 



11 

p+1 



et on effectue le changement de variable u — x x ; . 

e 

On arrive a I'estimation 



Or 



{x,e,ea,y)^xRo + e^-^^+y^^+y + 



X- 



4 



e dak ' ' ' V dak d&k dctk J d&k ' 
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dQ 

et la derivee partielle est done equivalente a — — ~ ex'^Ql quand e puis x tendent 

dak 

vers 0. D'ou la nouvelle approximation 

e'^'h' -^(p+l)//(0) e Ql je--'^' {ue' f ^^^/^) du 



aussi equivalente a 



Notons Ik = J^,^ e u'^du, qui est une integrale bien definie. Pour passer du 
chemin 71 a 7;, au moins dans un voisinage de 0, il sulfit d'effectuer une rotation 
d'angle ^^^^3^. Done, en posant 

/ 2(1 -l)n 
u — V exp I i 



p+l 



on obtient que 



I e-^^^\^du = j e-^^^^exp (^^^f^) exp {f-^^) dv 

/ 2(A; + !)(/- l)7r 
= exp I— — 

V P + l 
D'ou le resultat final : 

2lm 



[ exp ( - [ t^f(t)dt + r(x)] ^^(x,e,ea,y)dx CkS^'^'^^e exp( (k + V 



p+l 



la valeur de Ck etant la constante non nulle 



/ 2?7r 
Ck^C',*Ikexpl-{k + l)^ 



3.3.4 Existence de solutions y(x,e) exceptionnelles qui restent bornees 
dans un voisinage de taille fixe du point 

Theoreme 18 Soit V equation 

ey' = {x^fix) + eg{x, e))y + h{x, e) + ey'^P{x, e, ey) + Q{x, e, ea, y) (ill') 

oil toutes les fonctions qui interviennent sont holomorphes dans un voisinage Dg 
de X — (avec /(O) 7^ Oj, holomorphes en £ dans des secteurs centres en 0, P 
etant aussi holomorphe en ey — 0. Q se decompose comme en (3.16) 

Q{x, e, a, ey) — a{x, e)Ro{x) + eRi{x, e, a) + eyR2{x, e, a) 

+ ey'^R^{x, £, a, ey) + x^Ri{x, e, a) , 
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avec 

p-i 

a{x, e) = ^ ak{e)x'' = ao{x) + ea{x, e) . 

k=0 

Dans ce cas, il existe une unique solution formelle (y — yn(x)e'^, a — 
^ fln^") telle que ?/n(0) est bien definie pour tout n. 

On construit p+1 domaines Di C Dq contenant chacun une des montagnes 
(pour le relief 3?(/^ ^^^j—dt) ) naissant au voisinage du point col 0, chacun des 
domaines Di etant accessible a partir de son sommet, qu'on pourra eventuellement 
choisir infini (cf. fig 9). 

Alors, si pour un certain 6 fixe et tout e assez petit 

\9M\ + \h{x,e)\ + \P{x,e,eS)\ ^^^^ ^^^^ 

et 

^^"^'^l^^ est bornee sur (j D^, (H2) 
X j{x) 



1=1 



il existe, pour tout s > assez petit, des ak{s) et une fonction y{x, e) holornorphe 
bornee sur (J 2); et sur un voisinage complet de x — 0, B{0, p), tels que V equation 
(ill') soit verifies. 

Cette solution (y.a) correspond a la solution formelle, c'est-a-dire qu'elle ad- 
met les series formelles solutions comme developpements asymptotiques : 

oo oo 

y{x, e) yn{x)e'', d{£) ~ ^ a^e'^ . 

n=0 n=0 

Pour tout a dans un voisinage de taille e de ai{x) (trouve formellement, 
comme Oq precedemment), on sait que sur chaque domaine D;, il existe une 
solution e). Si les hypotheses du theoreme 18 sont verifiees, on a montre 

(au theoreme 14) qu'on pent prolonger ces solutions jusqu'en x = 0, et qu'elles 
restent bornees en ce point quand £ — 0. Pour demontrer le theoreme, nous al- 
lons montrer que pour tout e assez petit, il existe un a{x,£) tel que les fonctions 
?/['^(a;, £:) definies chacune sur leur montagne Di soient toutes egales en x = : 
?/W(0,£) = i/[2](0,e) = -.. = y[*^+il(0,5). 

Regardons comment y varie en fonction de a. L'equation (ill') implique, si on 
la derive par rapport a ccj : 

dy\' dQ . , dy 

x,e,ea,y) + —- x 



dak J dak dai 



xPf{x) + e [ g{x, e) + 2yP{x, e, sy) + ey "^(^^^ ^y) + ^ "^(^' ^' 
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Fig. 9: Exemple de voisinage de 0, avec representation de ©i et ©3 (p = 4) 



On voit qu'on a une equation differentielle lineaire en ^ ; les solutions peuvent 
done s'eerire : 

dyV] 1 



avee 



r{x) 



g{t, e) + 2yP{t, e, ey) + sy^—{t, s, ey) + -^(t, £, y) ] dt , 



dy 



e dy 



dQ 



une fonction bornee quand e tend vers 0, puisque est un O(e^). On pent 

dy 

appliquer le lemme 6, qui nous donne alors : 



dak £-»o 



p+1 



(3.22) 
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En reprenant le resultat complet du theoreme 17 (cf. (3.21)), on pent ecrire 
une formule de Taylor pour chacun des y^'l au voisinage de la fonction yo en a; = : 

t/W(0, £, 5) - t/o(0) = e%+m + e V afcC,£''+' exp l{k + 1)^ ] + o{. . .) 
qui donne, apres division par e, 

- (yM(0, £, ^) - yo(0)) = yp+i(0) + J] a.Cfe exp l^l{k + j + «(■ ■ ■ ) , 



ou ak = a/c^'*^^^. On rappelle que ni l^+i ni yo ne dependent de /, e on a. 
Nous allons regarder le vecteur 

2/P](0,£,fi)-|/o(0)) -r? 



1 f„[2]f 



[l (^y\P+^]{0,e,5^)-y,{0))-rj) 



de taillep+1, dependant desp + 2 variables e, wi —r],W2 — ckq, . . . , tWp+i = ctp-i. 
D'apres ce qui precede, la matrice jacobienne de y par rapport auxp+1 dernieres 
variables est 



dy_ 

Q^e=0,w={%+i{0),0,...,0) 



-1 Coe^^-^-?¥Tj 
-1 Coe(^'^'l?T) 



-1 Coe(^+^>^-|?^) ... C,_ie((^+^>^-|¥^). 



Le determinant correspondant est non nul. En effet, en factorisant le jacobien 
obtenu par colonnes, on obtient 



k=0 



exp 



exp 



2i7r 
p+ 1 



2m V 



1 I exp 



p+1 

2{p + l)i7r 
p+ 1 



exp 



exp 



p + 1 



exp 



p+1 

2i{p+l)n^^ 
p+1 



Le determinant qui reste est un determinant de Van der Monde, qui est bien non 
nul. 
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On en deduit que y est localement inversible autour de w = 0. 

On regarde la fonction y{e, w) ; on sait que, avec wq — (y^+i(0), 0, . . . , 0), 

y{o, wo) = 

et 

dy 



dw 



{e=o,w=wo) est inversible 



D'apres le theoreme des fonctions implicites, il existe done une unique fonction 
w{e) telle que y{e,w{e)) = 0, pour tout e dans un voisinage suffisamment petit 
de 0. Les fonctions a{e) que Ton en deduit sont bornees en e ; done les fonctions 
ea = ae'^~''~^^ tendent vers avec s. Cette propriete nous sufiit pour pouvoir 
effectivement appliquer le theoreme des fonctions implicites : les ckk sont toujours 
accompagnes d'un e en facteur dans toute ce paragraphe. Les fonctions ea{e) 
sont holomorphes en e dans des secteurs ouverts en 0, puisque le passage de £ a 
e' = ep+T n'est pas bijectif dans un visinage du point 0. 

Avec a, on trouve en meme temps la fonction y{x,s) correspondante a partir 
de y{0,e) = yo{0) + erj{e) ; dans les paragraphes precedents, on a montre qu'il 
s'agissait d'une fonction bornee, en a; = et dans tons les domaines autour de 
0. On vient ainsi de finir de demontrer I'existence de la vraie solution {y,a) du 
theoreme, bornee quand e tend vers 0, dans I'union des domaines 2)^. 

II reste a justifier que la solution formelle {y, a) est le developpement asympto- 
tique de {y, a). La demonstration se fait de maniere analogue a celle du theoreme 11. 
On sait que la solution {y, oq + sa) existe, avec ea bornee en e, et on ecrit 

y{x,£) = yo{x) +ez{x,e) 
ao{x) + ea{x, e) — ao{x) + eai{x) + s'^a{x, e) . 

Si on introduit ces notations dans I'equation differentielle (ill'), cela donne une 
equation differentielle lineaire en z avec un parametre a{x, e) qui intervient lui 
aussi de maniere lineaire : 



ez' = (x^fix) + eg(x, s)) z + ea(x, e)Ro(x) + x^ — — ^ ^ + 

e 

^ h{x,e)-h{x,0) _ ^ ^^^^ ^^^^ ^ ^^p^^^ ^ ^^^^^ ^^^^ ^ ^^^^ 

— * 

(on rappelle que itio ne depend que de x ; R\, R2 et R4 de x, e et eoi\ et que R^, 
depend en plus de ey). 

On verifie tres facilement que les hypotheses du theoreme sont encore verifiees 
pour cette equation : 
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1. f et g ne changent pas. 

2. la nouvelle fonction P est nulle. 

w — * 

3. la nouvelle fonction Q{x,e,ea,ey), quant a elle, se reduit a la somme des 

^ RAx,e,ea) — RJxjOjO) 
deux termes ea[x,e)Ro[x) + x'^ . Le premier terme 

ne pose aucune difRculte. Et le second peut s'ecrire aussi x'^R/i{x,e,ea). 

4. h{x, e) correspond a la deuxieme ligne dans I'equation ci-dessus. Cette fonc- 
tion qui peut s'exprimer en fonction de a; et £ uniquement (puisque y et ea 
sont bien connus) verifie la condition de majoration voulue. 

Avec ce qui vient d'etre demontre, on sait done qu'il existe z{x,e) et ea{x,e) 
bornees en e telles que z existe et est holomorphe en x — 0, et ea{e) est holomor- 
phe en e. 

Une simple recurrence suffit ensuite pour clore la demonstration du theoreme 18. 
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Chapitre 4 

L 'equation du Brusselator 



Le systeme differentiel du Brusselator apparait pour I'etude de revolution 
cinetique de certaines reactions chimiques auto-catalytiques. II pent se mettre 
sous la forme suivante : 

( ex = a{l + x)^y + a-x^ + (a - l)x 

1 y = -x(i + x) -2/(i + x)^ 

Qualitativement, il y a des solutions ncanardz de ce systeme pour des valeurs du 
parametre a voisines de ao = 1. Dans ce cas, la courbe lente a pour equation 

2 

y = (r+x)^ ' trajectoires des solutions peuvent etre amenees a suivre cette 

courbe lente y compris sur sa partie instable (cf. figure 10). On s'interessera en 
fait aux grands canards, c'est-a-dire aux trajectoires qui suivent la partie instable 
de la courbe jusqu'en x infini. 



4.1 Existence de solutions canards 



Nous allons commencer par transformer le systeme pour obtenir une equation 
sous une forme «normale», correspondant a celle etudiee dans la partie prece- 
dente (en particulier, nous ramenerons le point col pour le relief situe a I'infini 
en un point fini). Puis nous verrons quel est le domaine d'existence des solutions 
surstables a I'equation obtenue. 

Pour etudier ces solutions, nous allons prendre la nouvelle variable z — x, ce 
qui ramene la premiere equation du systeme a 

a(l + xfy = -«X^ - (a - l)x + ez 



et la seconde a : 

dy dy 2 , « - 1 £ -X - 

ax at a a a 
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-4 



-5 



-6 



-7 



-10 



-5 







10 



15 



Fig. 10: Courbe canard du Brusselator 



Or, la premiere equation ci-dessus nous donne aussi 



-2x a-1 1-x ^l2z + (l + x)dz/dx 



dy_ ^ 

dx (l + x)3 a + ' a (l + x)2 
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soit 

a{l + X) -J- ^ -a + r;—, 'rSl— h 3- . 

dx 1 + x \i + x dxy 

Le systeme du Brusselator peut done se ramener a I'equation suivante, que nous 
appellerons equation du Brusselator. 

sz— — z — - — \ + [a — l)z + z \ z — ' 



dx 1 + xV J \ 2 Jl + x 

On voit deja sur cette equation que la fonction $o(x) = ''^"'2^'' est particuliere 
pour cette equation : pour le parametre a = 1 elle est solution quand e = 0, et 
c'est alors la seule solution (pour £ = quand on fait varier a) qui soit continue 
en X = 0. 

On effectue le changement de variable suivant, qui ramene I'infini en —1 et 
laisse en 0. 

X = ^ — , puis — = -(1 + x) 

1 + x dx 

Si on compare alors avec I'equation de Van der Pol, correspondra ici a 1 (le 
point col qui n'est, exceptionnellement, pas un point tournant), —1 a —1 (I'autre 
point col du relief). 
D'ou I'equation 



dz{x) 2x / 1 \ a — 1 / 1\ 



2s 



(l + x)n 2(l + xfJ (l + x)2 V 2(l + xfJl + x 

(4.1) 

Montrons que cette equation est bien du genre etudie lors de la generalisation 
(voir partie precedente). On pose pour cela z — ^o{x){l+ey), ou $o(a^) correspond 
a la courbe lente, ici 2(14^5 remplace dans I'equation (4.1) apres en avoir 

divise les deux membres par z 

-2 2x61/ a~l ^ 2e 



£$'o(l + ey) + eY<l>o = - - $0^2/: 



[l + xfil + ey) (l + x)2 " ^1 + x 
Done 

e 



4x(l + x) + 



1 + x 



H 2{x + l)- — - -ey'^4:x{l + x)- ^ 



1 + x e 1 + £y 



On voit que autour de 0, on a : 
• p = 1 et f{x) = 4(1 + x) 
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• h{x, s) — 

• P{x,e,ey) 



3 
1+x 



4x{l+x) 



• Q{x, e, y) = —2{x + 1) avec a — 

• On prend alors pour T)q I'ensemble C prive du demi-axe reel ] — oo, —1]. 

• Toutes les hypotheses etant respectees, vu le rehef represente figure 11, le 
theoreme 18 s'applique : pour une valeur de a bien choisie, soit a'^{s), il 
existe une solution holomorphe z~^{x) qui sera continue en et qui a pour 
domaine d'existence les montagnes Nord et Est et les vallees adjacentes ; voir 
fig. 11. Symetriquement, en prenant les conjugues, on trouve une solution 
z~ — z+ correspondant au parametre a — a~{e) — a+. 

Bien evidemment, ces deux solutions (a, z) correspondent toutes les deux aux 
series formelles uniques 



construites de telle sorte que tons les termes Zn sont continus en a; = 0. 

Si on regarde ensuite ce qui se passe autour de —1, en changeant de variable : 
X — X — 1, on trouve de maniere analogue 



• On prend alors pour D la montagne Est (a partir de —1) et eventuellement 
une partie des vallees adjacentes. 

• Les hypotheses des theoremes 12, 13, 15 sont respectees : on pent prolonger 
le domaine jusqu'en — 1 + Xi^fe G M ; et la solution canard admet une sorte 
de developpement asymptotique en ^/e pour des x tels que |l + a;| < Iv^l^"*"^ 
{y > 0), obtenu a partir de la solution formelle de I'equation difi^erentielle 
non singulierement perturbee ; tandis que le developpement asymptotique 
z{e. x) = 2;„(,t)£" reste valable, dans le sens du theoreme 15, pour les x 
(sur la montagne Est) tels que |1 + a;| > < u < /j, < 1. 

4.2 Estimation de a+ — a~ 

L'existence de solutions la ou on le souhaite est done acquise. Maintenant, 
pour des raisons de simplicites de calcul, nous allons repartir de I'equation (4.1) ; la 
progression suivra alors exactement le meme shema que pour I'equation de Van der 
Pol. C'est a dire que nous allons etablir une equation reliant la difference — a" 




• p = 1, f{x) = 4(1 - x) et e' ^ ^/e 

• g{x, s) — ^ d'ou aussi Gq{X) — 

• h(x,s) = I - 2x^ d'ou aussi Ho{X) = 
. Pix,e,ey) = d'ou aussi Po{X,Y) 



3_ 
X 



4X 
1+Y 
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Fig. 11: Relief correspondant a 2t{l + t)dt : domaine accessible a partir des sommets 
Est ou Nord, pour a = a'^ ; les courbes passant par sont les courbes de niveau R{x) = 0, 
celles passant par —1 ont pour equation R{x) = 1/3. 




a la difference entre les fonctions z'^ — z~ , etudiee au voisinage du point singulier 
— 1 ; pour pouvoir ensuite estimer chacun des termes presents dans I'equation 
obtenue, nous aurons besoin en particulier des equivalents (pour e tendant vers 
0) de toutes les fonctions solutions des equations interieures comme exterieures, 
et des coefficients de Stokes des solutions interieures. Certains des calculs seront 
fails ici en preliminaires. 

Pour etudier le voisinage de —1, le changement de variable qui s'impose est en 
tons cas le suivant 

D'ou I'equation interieure, non singulierement perturbee 

_ 2(-l + ViX) / 1 X g-l 21-/ 1 X 
dX X2 V 2X^) X^ X\ 2X^J ■ 

Cette derniere equation a pour solution formelle la serie Y — ^y„(X)£"/^ (et 
on sail que cette serie correspond a une vraie solution dans certains domaines), 
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avec lo solution de I'equation obtenue en posant £ = 0, soit 

4.2.1 Quelques equivalents 

Si on cherche des equivalents pour les premieres fonctions des series y et F, 
on commence par ecrire dans I'equation (4.1) la fonction z et le parametre a sous 
la forme z{x) = ^o{x) + ^^i{x) + o{e) et a = 1 + aiS + o{s) ; on part de 

fl + y)^ 1 



2(l + a;)3 ' 
puis, comme 

2y 
1 + X 

et que $i doit etre continue en 0, on trouve 

ai = 3/2 , $i(x)- ^ 3 ^ (2 + x) 

autrement dit 

^ t ^ 3 2 + a; 
$1 X = 



8(l + a;)5 ■ 

On continue pour obtenir le terme suivant : dans I'equation du Brusselator, les 
termes de degre 2 en £ meneront a I'equation 

1+X 1+X 
d'ou ^ 

02 = 15/8 , Mx) - ^^^^/^ (6x' + 29x' + 48x + 30) 
autrement dit 

^ , , 3 5a;3 + 23^^ + 42a; + 30 
^^(^^^32 (iTW ■ 

Ensuite, pour la serie F, on fait le changement de variable et on identifie, avec 
X = —1 + Xy/e, les deux series 

e^/^[^o{x) + £$i(x) + e^^2{x)') = Yo{X) + VEY,{X) + eY^iX) + o{e) , 

ce qui mene a 



X-.0O 2X3 gX5 16X7 ' / gjs^4 32X6 ' ^ 4^5 
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4.2.2 Coefficient de Stokes pour I'equation (4.2) 

Si I'equation (4.2) a deux solutions conjuguees Yq' et Y^ , bornees toutes deux 
pour les X e R assez grand, elles out necessairement la meme serie asymptotique 
(dont les termes sont donnes ci-dessus) et la difference entre les deux est done 
exponentionnellement petite. Elle correspond a des termes de Stokes ; nous allons 
montrer que : 

Lemme 19 La difference entre deux solutions bornees en I'infini de I'equation (4.2) 
est equivalente, quand X tend vers I'infini, d 

(Y+ - 17) (^) ~ 32iv^X^e-^^' . 

Pour parvenir a ce resultat, nous allons essayer de donner, sous une certaine forme 
(car on ne pent esperer en donner explicitement les solutions exactes), les solutions 
de I'equation. 

(4.2)yo' = -4^o + 4i-4+ ^ 



On commence par changer de variable en posant t = 1/X 

dYn '^Yn 9 t"^ 

dt ~ t ^ Fo' 
puis on introduit la nouvelle fonction u{t) telle que lo(0 — t^u{t) 

2tu + t^u' =2tu -t^ + 2-t/u 

et v{t) verifiant u{t) ^ -t - 2r ^ + v{t) 

dv 1 



dt t'^ + 2- tv{t) ■ 

A present plutot que de regarder v en fonction de t, on prend la fonction re- 
ciproque : on cherche t en fonction de v 

^ ^t{vf + 2 - vt{v) 
dv 

On reconnait une equation de Ricatti. La maniere habituelle de resoudre ce type 
d'equation est de poser t — 

z z ^ z z 

+ ^^2 + -v + — 

z z^ z z^ 

z" + vz' + 2z =0 
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Cette derniere equation est une equation differentielle hypergeometrique, qui a 
les deux solutions independantes suivantes (parmi d'autres possibles) : 

de derivees 

z[(v)^(l-v^)e-''"/^ et z'Jv)^(l-v^)e-''"/^ + - . 



On pent finalement ecrire toute solution t{v) sous la forme : 



^ ^ J too ■w' 



Mais toutes ces solutions t ne nous conviennent pas. En effet, on a vu precedem- 
ment que les solutions Y considerees devaient verifier la propriete suivante : 



X-.00 2X^ ' 
soit, apres changements de variables, 

-t^ -2t + fv - , 
ce qui suppose au moins que 

t{v) - 2/v. 

v—*<x> 



Or si on prend C2 = 0, on trouve t — qui ne verifie pas cet equivalent. 

Par contre, pour toute autre solution (dans le secteur ou elle est definie a 
I'infini), on trouve le bon equivalent. On montre en utilisant des integrations par 
parties que, quand v ^ 00, 

Jioo W \ ) 



Ce qui fait que, si C2 7^ 0, 



A^) ~ (4.3) 
z{v) ~ \ (4.4) 

v—*oo V 
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et done 



2 

t{v) ~ - 

v—*oo V 

ce qui donne une solution Yq qui se comporte bien comme prevu a I'infini. 

Comme C2 7^ dans les cas qui nous interessent, on ecrira desormais t{v) avec 
une seule constante C — C1/C2 : 

t{v) = ^ ^ 



ve 



\Jtoo w' I 



Si alors C = 0, la fonction t~^{v) est bien definie pour tout v tel que aig{v) e 
[— 7r/4 + 6, 57r/4 — 6],S > 0, parce que pour arg(u^) g] — 7r/2, 7r/2[, I'exponentielle 
e"" est bien bornee et que pour arg(i') e [7r/4, 37r/4], ce sont les integrales qui 
sont effectivement exponentionnellement petites. Si on repasse a la variable X, 
cela donne I'existence d'une solution au moins pour tout X assez grand sur 
I'axe reel et la montagne Nord (cf. figure 11) d'apres I'equivalent X = l/t{v) ~ 2v. 

Pour C = iV27r, on pent reecrire t{v) sous la forme suivante : 



^— V / J —too V 



J —too w' 



En effet, en integrant par parties 

/ w w I ■ 

J —100 ^ L J —too —too 

et avec le changement de variable t — iw/ -\/2 , 

/•— 00 

= V2/i / exp{-t'^)dt = iV27r 
J +00 

La solution t~ donne une solution pour des X dans un domaine symetrique 
(par la conjugaison complexe) au domaine image de t~^. 

On considere a present des variables complexes v'^ et v' telles que 

t+{v+) ^r{v-) e R. 

Pour des raisons de symetrie, ces deux nombres sont a priori conjugues. On 
souhaite etudier la difference entre ces deux nombres quand ils tendent vers +00. 
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La methode sera semblable a celle utilisee pour I'equation de Van der Pol. 

t+(v+) =t-{v-) 
z~^'{v~^) z~'{v~) 



z+ (v-) il+i 



Comme les termes e^'"^ 1'^ jz^^ sont exponentionnellement petits quand v — > cxd, on 
pent ecrire des developpements limites : 



z' \v 



z^[v~ 



1 + i^27r 



t+{v+) - t-{v-) r.t+{v-)iV2^e-''"/^ 



(^;2 _ i)e--V2 ^;e-V2' 

z~^'{v~) z~^{v~) 

{v^ - 1) V 
z+'{v-) ~ z+{v-) 



+ o{...) 



D'ou, avec les equivalents pour z et z' trouves en (4.4) et (4.3) 

t+(v+) - r(v-) ~ - i^e-^'/V 

Or, comme v'^ et v~ sont proches, on pent aussi ecrire 

t+{v+) - t-{v-) - {v+ - V-) t+'{v+) . 

Done finalement 




On en deduit facilement les equivalents correspondants pour {Y(^ — Yq ){X), 
sachant que 

v+ -v~ 2 
yo^-yo-=^^ et v^-^2X. 

On obtient 

{¥+ - Yo-){X) ~32^y2^XV2^' , (19) 
avec done un coefficient de Stokes 

C = 32iV27r. 
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4.2.3 L'equation differentielle pour I'estimation de a+ — a 

Comme dans la premiere partie on cherche une equation differentielle en d(x) 
(2;+ — z~){x) ; on pose aussi b = a'^ — a~ . 
A partir de l'equation (4.1) ci-dessous 



dz+{x) _ 2a; / 1 \ a+ - 1 / + 1 \ 2£ 

^ dx ~ (l + xY\ ~ 2z+(1 + xY) ~ + ~ V ~ 2(1 + xY)T 



{i + xy\ 2z+{i + xYJ {i + xY V 2{i + xyJi + x 

et de la meme equation ecrite pour {z~,a~), on arrive, en faisant la difference, a 



{l + x)^\z- z+J {1 + xy 1 + x 
, fx 2e \ b 

ed'{x) = d (^(^^^^5^-^+ Y+^J ~ {l + xf ■ 

D'ou, pour tout xq reel, 

d{xo) = — / 



b n exp^-— 

^-dt 



avec G't(x) = (i+t)5z+(t)z-{t) '^^ ~ 2e\n{l + x), soit 

4.2.4 Calcul de I'equivalent 

On peut estimer chacun des termes du produit de (4.5) (on travaillera ici avec 
un e reel positif). 

Avec la methode du point col. 



J +00 



' +00 

Puis 



g-G(t)/£^^ V2^\z{0)\ avec \z{0)\ ~ $o(0) = J • (4.6) 



d{x^)^e-^^l^{Y^{X,)-Y^[X{)). (4.7) 

Reste e'^'^^")/'^, ou on prend .tq = — 1 + Xi^. On coupe I'integrale corre- 
spondant a G{xq) en deux, en un point dependant peu de comme par 
exemple \[e) = —1 + ; le but etant qu'a droite on puisse utiliser I'estimation 
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z{t,e) ^^o{t)+e^i{t)+e'^^2{t) + o{e'^-'^/^^) (theoreme 15) ; et a gauche, grace au 
corollairelS, restimation s) = -{Yo{t/e)+^Yi{t/e)+eY2{t/e))+o{^^~^^^°). 



e 

••xo 



f t t 

Pour Xi = (i+ty'z+{t)z-(t) ^'^ utilise done pour z la forme z = $o + + 
o{e) (on a bien o{e) car $2 est bornee), qui donne 



z+{t)z-{t) 



done 



= + 2A^ - £(3(1 + A)=^ - 3 - 61n(l + A)) . {KIi) 

Pour 22 = J_^° {i+tfz+{t)z-{t) place au voisinage de i = —1, et on 

change de variable : t — —1 -\- X^Je, z — e'^/'^Y . 



A+i Vi'^%-3r-r+ x^Y-{x)Y+{x)' 



ou on pose bien sur £ = Le calcul ressemble a celui effectue plus haut dans 
le cas de van der Pol ; on commence par developper les fonctions dans I'integrale 
suivant £, avant d'estimer chacun des termes obtenus, en verifiant que le dernier 
terme, et les suivants, sont bien negligeables par rapport aux precedents. 
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Comme Y(X) ~ 2x3^ I'integrale converge, y compris pour Xq — — 1, et 

e r -^ + v^-^ dX - e'l- r -^ + ^^ , dX 

% X^Y,-{X)Y,^{xf I X^Y,-{X)Y,^X)W Yo-J 

^ V ' ^ V ' ^ V ' 

(A) (B) (C) 



+ 2e' 



(D) 

On s'occupe de chacun des termes de la somme 

—o{e) ( Rappel : £ — — ^ ) 
V V £ / 

-£(-6(A + l) + 3(A+l)') 



, ^3/2 ^ ^4^2 _ g ^ ^ (^/^^^ ^ _^3/2 _ 



6£ 



_^(A + l)3 + 6£(A+l) 



Pour A, il faut etre plus precis et distinguer Y^ de 1^ ; on repasse done par 
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I'equation differentielle (4.2), 



2 .... .... . 1 -V + V 



(V-n-)' = -v(n^-n-) + 



qui donne, en integrant I'equation par la methode de variation de la constante, 
pour tons ^, positifs : 



D'apres le resultat (19) obtenu plus haut, 

A ^ £ (in - Y,-) (Xi) - In (cfe-''^) + 2 In ^) 

-e^^ln {Y+ - Fq') (Xi) - \nC - 41n(A + 1) + 2^^^^ + 21n£ 

+ 21n(xo + l) -21n(A + l) ) . 



Enfin, 



I2 ~ A+B+C+D ~ -^A' - 2A' + ^ + £ (^In {Y^-Y^-) {Xi) + 2 ln(a;o + 1) 

-61n(A + l) + 21n£-lnC + 3(A + l)M , (i^Xs) 



puis en reprenant le resultat {RIi) 

exp (^^) - exp (1) (ro+ - (X,)(xo + if^e' , 
done avec (4.7) et (4.6), I'equation (4.5) devient 

-3/2.^+ y-VY^ , V^e3 (ro+-yo-)(x,)(a;o + i)W/3^ 
s [Y, -n ) (XO ~ 6 2CiKTTF 

3/2 ^ ^eV2^^3/2g^ 
2C 

^^^cxp(-4) 
£3 
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Finalement, on arrive au resultat du 
Theoreme 20 

(a"*" — a~){e) ~ 64ie~^exp (^~'^ ) ^' 

4.3 Consequence pour le developpement en serie de a{£) 

Corollaire 21 Les coefficients an du developpement asymptotique de a+ et a~ 
ont une croissance Gevrey, et verifient, quand n ^ oo, 

an ~ 54n' n! 

Etudions d'abord ce qui se passe pour la fonction a~^{e) quand on fait varier 
I'argument 9 de e. II est clair que a'^{s) et v^{x, s) dependent holomorphiquement 
de e dans des secteurs centres en (mais pas en lui-meme). Regardons en 
particulier ce qui se passe quand on fait varier 9 entre a — 27r. En etudiant le 
relief correspondant a 9, on pent deduire le domaine d'existence de la solution 

correspondante (qui varie continument avec 9 lui aussi) ; cette evolution du 
domaine est represente figure 12. On remarque tout de suite que pour e reel, le 
domaine de pour ee~^^'^ correspond au domaine de z~ pour la valeur s. Par 
unicite des solutions, on en deduit que a"'"(£e~^*'^) = a~{e). C'est de cette propriete 
que Ton se servira ci-dessous. 

On pent faire exactement les memes deductions que pour I'equation de Van 
der Pol, les deux solutions y'^{x,e) et y~{x,e) etant reliees de la meme maniere. 
On integre a nouveau sur un chemin autour de £ = (celui de la figure 6), pour 
trouver que 

^"^2^70 
En utilisant le theoreme 20, 

~ — /" ° GAie-^exp ( -—] £-V"-^(i£ 

32e-3 n f 2\ 

exp -— ] e " de. 



Puis on pose i = 2/3£ et on transforme I'integrale en une autre, equivalente : 



32e-3 /•+°° /3^ 



/ \ 9 

32e-3 /3\"+3 



e-H^'+^dt 



anr^—^i-) (n + 2)! 
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Fig. 12: Domaine d'existence de z'^ pour differentes valeurs 6 = arg(e). 




Ill 
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